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Abstract 

Symbolic codes for rotational orbits and "islands-around-islands" are constructed for the 
quadratic, area-preserving Henon map. The codes are based upon continuation from an anti- 
integrable limit, or alternatively from the horseshoe. Given any sequence of rotation numbers 
we obtain symbolic sequences for the corresponding elliptic and hyperbolic rotational orbits. 
These are shown to be consistent with numerical evidence. The resulting symbolic partition 
of the phase space consists of wedges constructed from images of the symmetry lines of the 
map. 



1 Introduction 

Symbolic dynamics has been profitably used in the study of many dynamical systems since its 
invention by Hadamard in 1898 and naming by Morse and Hedlund in 1938 0121 El- Among its 
many uses, symbolic dynamics can provide useful information about topological invariants such 
as the enumeration of periodic orbits and the entropy. It can also lead to numerical methods 
for finding periodic, homoclinic and chaotic orbits, and facilitate characterization of transport. 
Symbolic dynamics applies most directly to systems that are hyperbolic; indeed it was invented 
to describe geodesies on surfaces with negative curvature and famously applies to hyperbolic 
toral automorphisms and the Smale horseshoe. 

In this paper we continue the study, began in jlj, of the coding of orbits of Henon's quadratic, 
area-preserving mapping [SI. Previously we used the concept of an anti-integrable (AI) limit [S] 
to define codes by continuation (see 321) and studied the bifurcations of homoclinic orbits that 
destroy the horseshoe of this map. In the current paper we study the codes of orbits that are 
born in rotational bifurcations of the elliptic fixed point. The rule that we obtain identifies the 
subset of orbits in the horseshoe that become rotational orbits, encircling the elliptic fixed point 
of the Henon map. 

The extension of symbolic dynamics to systems with stable orbits has proved difficult, ex- 
cept for the case of one-dimensional maps where classical results apply |3 |H1 • The dissipative 
Henon mapping 9^ is a natural system to attempt the generalization of these results to multidi- 
mensional, nonhyperbolic systems, especially as it reduces (when 6 = 0) to the one-dimensional 
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logistic map. Just as the symbols for the logistic map are based on the itinerary of an orbit rel- 
ative to the critical point, Grassberger and Kantz [TOj proposed that the symbols for the Henon 
map could be determined by partitioning the plane using the stable and unstable foliation and 
their "primary homoclinic tangencies," that is, points at which the local curvature of the un- 
stable manifold diverges. This led to the concept of a "primary pruning front" in the symbol 
plane as determining the allowed sequences for a particular map . Symbol sequences defined 
in this way can exhibit monodromy, that is, morph into new sequences on paths that encircle 
codimension-two bifurcations, such as a cusp along these paths, the primary tangencies 
may exhibit discontinuities jEl- As far as we know, this method has not been applied to the 
area-preserving case (6 = ±1), though it has been applied to other dissipative systems, such as 
the cubic, generalized Henon map |14j . 

Several other methods have also been proposed to obtain symbolic codes for orbits of the 
Henon map. Biham and Wenzel defined the codes as the signature of a pseudo-gradient method 
for finding periodic orbits ^Hl- Sterling and Meiss proved that this technique works sufficiently 
close to an anti-integrable limit jT^. Unfortunately, it does not always converge to fixed points, 
and sometimes gives two codes for the same orbit 17^. Hansen and Cvitanovic defined codes by 
approximating the two-dimensional map near 6 = by a sequence of one-dimensional unimodal 
maps |18) . The code defining an orbit can change if it crosses the critical point of one of the 
approximating maps. 

Primary homoclinic tangencies have also been used to construct codes for the area-preserving 
standard map. For parameter values where elliptic periodic orbits and their associated islands are 
small (large k), Christiansen and Politi have shown that a primary set of homoclinic tangencies 
can be identified by their proximity to the dominant fold lines in the map P^. Gaps between 
these points can be connected with symmetry lines associated with the reversibility of the map 
to form a curve that creates a symbol partition [2111 • The symbol boundary can be modified to 
include elliptic islands, |211 l2Uj . by choosing appropriate images of the symmetry lines. This 
method does not explain why symmetry lines are important, nor does it give a recipe for selecting 
the proper lines. 

In ^we give a rule for constructing the symbolic codes for orbits of the Henon map with any 
given rotation number. These codes are consistent with those obtained by numerical continuation 
from the AI limit. We also show that these codes have definite symmetry properties, and that 
the symbolic partition in an elliptic island has the form of a wedge with apex at the elliptic 
fixed point and whose boundaries are constructed from specific symmetry lines. These wedges 
are similar to those found by Christiansen and Politi for the standard map. 

Our rotational codes are closely related to those for maps with a natural angle variable, for 
example for circle maps |22| \2'A\ and cat maps |24j. We review these ideas in the Appendices. 

We also develop a systematic rule for obtaining the symbolic codes of rotational orbits with 
"higher class" |25j . that is to "islands-around-islands," in ^4151 These correspond, for example, 
to orbits that rotate around orbits that rotate around the elliptic fixed point; thus they are 
defined by a sequence of rotation numbers. Our construction provides codes for the elliptic 
and hyperbolic orbits for each sequence of rotation numbers. Again the codes are shown to be 
consistent with the numerical evidence. The resulting symbolic partition is constructed from a 
sequence of wedges defined by symmetry lines. 

A different method for constructing symbolic codes for islands-around-islands was given 
previously j26| I27j ; however in these cases the entire set of orbits in an island was assigned the 
same sequence and the motivation was to study the transport implications for chaotic orbits 
outside the islands [23 [2Hj • 
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2 Symbolic codes for the Henon map by continuation 



The area-preserving Henon diffeomorphism can be written as 

(x, x') ^ f{x, x') = (x', -x-k + . (1) 

An orbit of this map is written as a bi- infinite sequence O = {. . . x^i, xq, xi . . .) with (x^+i, = 
f{xt, The dynamics of this map exhibit the full range of behavior of typical area-preserving 

maps, including infinite cascades of period-doubling bifurcations, invariant circles, cantori, trans- 
port, twistless bifurcations, etc. [1^11301 1^- As k increases, the dynamics become increasingly 
chaotic, and beyond a critical value, k = kn, the set of bounded orbits forms a Smale horseshoe 
(this set has measure zero — almost all orbits escape to infinity). It was proved by Devaney and 
Nitecki that there is indeed a hyperbolic horseshoe when A: > 5 + 2-v/5 > kjj \3'Zl I16j . Our 
numerical studies indicate that ^1 11] 

kn ~ 5.699310786700 . (2) 

The dynamics of the bounded orbits in the hyperbolic horseshoe is conjugate to the full shift 
on two symbols. A simple way of constructing these symbols for the Henon map is to define 

st = sgnxt (3) 

because the bounded invariant set is divided by the coordinate axes, see Fig. ^ 

An alternative method for obtaining these codes is to view the Henon map as arising by 
continuation from an anti-integrable (AI) limit 01331; © define 

z = ex , with e = —== , (4) 
Vk + l 

to rescale the map as 

e{zt+i + zt-i - e) = zj - I . 

To define e, we assumed that A; > — 1, but this is not much of a restriction since there are no 
bounded orbits for the Henon map when k < —1. The case e — > is the AI limit [161 In 
this limit orbits reduce to sequences zt = ±1, where the choice of sign is arbitrary — the map 
is equivalent to the full shift on the symbols it. Every such sequence continues away from the 
limit to an orbit of the Henon map. Conversely, when there is a horseshoe these symbols agree 
with Eq. Q, and every bounded orbit of the Henon horseshoe continues to an orbit at the AI 
limit. 

Each point in the horseshoe is represented by a bi-infinite sequence of signs, together with 
the binary point representing the current position: 

{xt,Xt+l) = S = (. . . St-2 St-l.StSt+l . . .) 

We denote the right-shift by cj, so that 

as = {...st-iSfSt+i...) (5) 

When there is a hyperbolic horseshoe for k > kn, the Henon map /, restricted to the set of 
bounded orbits, is conjugate to the shift map a acting on the sequence s that is the code for x. 

^We do not use Henon's original form, (^,?7) ^ (1 — a^^ + r},bS^), since it becomes linear at a = 0, and is 
orientation-reversing for 6 > 0. Our form can be obtained by setting b = —1, and defining k = a, x = krj, and 
x' = -k^. 



3 



Figure 1: Stable (blue) and unstable (red) manifolds of the hyperbolic fixed point for the Henon map when 
k = 6. The closure of the intersection of these manifolds is the horseshoe of bounded orbits. 

We use continuation away from the AI limit to find orbits with a given code as k varies |16j . 
Such a "global coding" is complete if every smooth one-parameter family in the extended space 

X ]R (phase space times parameter) connects to the AI limit. Because of twistless bifurcations 
it is not enough to simply take the map parameter A; as a family parameter [311 . We do not 
know any counterexample to this "smooth no-bubble conjecture", and adopt it as our working 
hypothesis. 

We use the signs -|- and — to denote the symbols st- For example the hyperbolic fixed 
point, at xt = 1 -|- Vl + k, has code st = +; we denote the bi-infinite sequence for this orbit by 
s = (-|-)°°. In general a period-n orbit is given by the bi-infinite concatenation of a sequence of 
n symbols; we represent this with a superscript oo: (sq si . . . Sn~i)°°- We often denote repeated 
symbols with a superscript, thus {—+'^)°° is a period-8 orbit. The parity of a finite symbol 
sequence is defined to be the product of its symbols: 



Equivalently the parity is even or odd if the number of minus signs is as well. A parity can be 
assigned to any periodic orbit by computing the parity of its fundamental sequence; for example, 



Using the continuation method, we can follow orbits from the AI limit to visualize the relation 
between the symbolic codes and the positions of the corresponding points in phase space. If we 
distinguish points solely by a single symbol, sq, this gives a partition of the phase space into 
regions corresponding to the -|- and the — code. Using the conjugacy between the shift a and 



n— 1 




(6) 



t=o 



7r((- + ++r) = -l. 
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the Henon map on these orbits, we could reconstruct the fuh symbol sequence for each orbit by 
following the sequence of visits to the two elements of the partition. 




Figure 2: Points on the 99 periodic orbits of the Henon map with periods up to 10 for two different parameter 
values. In the left panel A; = 24 (e = 0.2) and in the right k = 5.5 (e ~ 0.392). The blue (red) curves are the initial 
lobes of the stable (unstable) manifolds for the hyperbolic fixed point, (+)°°, located in the upper right corner. 

In Fig. [21 a point is colored cyan if sq = — , and magenta if sq = +. We use the scaled 
coordinates {z,z'), defined in Eq. (jlj for the plots, so that as A: — > oo all bounded orbits 
converge to the four AI states {z,z') = (ibl,ibl). As k decreases the trajectories move away 
from these points as shown in the left frame of Fig. [3 For this value of e it can be proven 
that all of the trajectories reside within the union of four small squares with sides M^o ~ 0.568 
that are centered on the four AI states ^01 — actually we observe that they are contained in 
smaller rectangles bounded by segments of the stable and unstable manifolds of the hyperbolic 
fixed point. These segments intersect at points on the "type one" homoclinic orbits with symbol 
sequences +°° . — (— ) — \-°° and +°° . — (+) — \-°° j^. As e increases the bounded orbits move 
further from the AI states, and the stable and unstable manifolds of (+)°° approach their first 
tangency. Just below k = kn (right panel of Fig. [21) the two type-one homoclinic orbits have 
been destroyed in a homoclinic bifurcation and only a subset of symbol sequences are now 
allowed. Nevertheless, the symbol boundary is still very simple (at least up to this period), and 
is delineated by the near tangency between the first lobes of the stable and unstable manifolds. 

One phenomenon special to area-preserving maps is the existence of elliptic periodic orbits. 
For example, the second fixed point, (— )°° (at xt = 1 — \/l -1- k) becomes elliptic at /c = 3. More 
generally this orbit has positive residue for any k > —1. Recall that Greene's residue, r, is a 
convenient encoding for the linear stability of an orbit of an area-preserving map: given a point 
(x,x') on a period-n orbit of /, then 

r^i(2-Tr(Z)r(x,x')) (7) 

Orbits are elliptic when < r < 1, hyperbolic when r < and reflection hyperbolic when 
r > 1. For example, the fixed point (— )°° is elliptic for —1 < A; < 3 and at /c = 3 undergoes a 
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period-doubling bifurcation becoming hyperbolic with reflection. More generally [3] for a period 
n orbit, as A; ^ oo, 

r~-vr(s)i(4A;r/2. (8) 

Thus the parity of a symbol sequence determines the character of the orbit for large k. Often, 
we refer to orbits that are born with positive residue as "elliptic" even though they typically 
undergo bifurcations that destroy their stability and may even smoothly continue to orbits with 
negative residue. 

Our goal for this paper is to investigate the structure of the symbol partition when there 
are elliptic periodic orbits and their associated rotational orbits. In particular. Fig. |21 shows 
that the symbol boundary near an elliptic orbit exhibits a structure radically different from the 
horseshoe shown in Fig. |2l In the following sections, we will explain the wedge shaped symbol 
boundary that eminates from the elliptic fixed point in Fig. IHI To do this, we must construct 
the codes for rotational orbits. 




-1 1 

z 



Figure 3: Periodic orbits of the Henon map for k — 2.0. Plotted are points on all periodic orbits up to period 
10 that still exist, as well as all rotational (class-one) orbits with periods up to 31 that exist at this parameter 
value. Points with so = + ( — ) are colored magenta (cyan). It is apparent that the first few lobes of the stable 
and unstable manifolds no longer delineate the symbol boundary in the vicinity of the elliptic fixed point. 



3 Rotational (Class-One) Orbits 

One special class of orbits for the Henon map are the rotational orbits — those that are created 
by a rotational bifurcation of the elliptic fixed point (— )°^. Generally, a rotational bifurcation 
occurs when the linearized rotation number at an elliptic fixed point passes through a given 
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rotation number < w < |, 



i.e., when 



r = 




) = sin^(7rc<j) 



(9) 



Since the residue of the elhptic fixed point is r{{—)°°) = ^ 
uj to obtain 



+ we can solve Eq. @ for k or 



k{io) = cos{2ttu))[cos{2ttu)) — 2) 




(10) 



Since the residue of a period-n orbits grows like k^/'^ as k 
numerical work to use the scaled residue 



(X) it is more convenient for 



1 



sgn(r) [4|r|]^/' 



(11) 



P = 



^/k+l 



Notice p = <^ r = 0, and as — > oo, Eq. ^ implies that p ±2. 

When uj{k) is rational, a rotational bifurcation typically^ creates a pair of orbits, one elliptic 
and one hyperbolic, we call these class-one rotational orbits following Meiss H5] (the two fixed 
points are therefore class-zero orbits). A theorem of Franks implies that there are class-one 
period q orbits for each rational rotation number p/q whenever < p/q < uj{k) [SHI- Each 
elliptic class-one periodic orbit is typically encircled by rotational orbits of class-two, giving rise 
to "island chains." We will discuss these higher class orbits in ^ When uj{k) is irrational and 
Diophantine, Moser's twist theorem implies that the rotational bifurcation gives birth to an 
invariant circle j^n]- If uj{k) is not sufficiently irrational, we expect that a cantorus is born (as 
indicated by Aubry-Mather theory, though we know of no proof of this for the Henon map). 

We begin this section with some observations from our numerical experiments. We then 
construct a rule that generalizes these observations and appears to be consistent with all class- 
one orbits. This construction has a natural geometric interpretation in terms of circle maps. 
The code that we obtain is related to that for rotational orbits in other systems, for example 
the codes obtained for minimizing rotational orbits in twist maps by Aubry for cat maps 
by Percival and Vivaldi [25) billiards by Backer and DuUin |3H], and for circle maps by Lin [2]. 
We discuss some of these relationships in the Appendices. 

3.1 Preliminary Numerical Observations 

It was not obvious to us which of the possible period-n orbits at the AI limit should be identified 
as the class-one rotational orbits. Thus we began by searching for periodic codes that continue to 
orbits that collide with (— )°° at k(uj) defined by Eq. A brute force way to systematically 

search for rotational orbits, would be to choose a period and continue all periodic codes of 
that length from the AI limit until the corresponding orbits bifurcate. For example, at period 

4 there are three distinct periodic codes: ( 1")°°)( 1"+)°°) and ( — 1-++)°°. Upon 

continuation, we discover that the orbits corresponding to ( 1-+)°° and ( — h ++)°° collide 

when k = k{j) = at the position of the elliptic fixed point. These two orbits have negative 
and positive residues, respectively (in agreement with Eq. Q) during their entire lifetimes, and 
thus they represent the hyperbolic and elliptic | orbits. In our previous numerical experiments, 

^When the period is larger than four. The tripling and quadrupling bifurcations are special, see e.g. |.34| . 
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we followed each of the 1.47 million orbits of the Henon map with periods up to 24 in order to 
obtain a numerical estimate of the parameter at which the Henon map horseshoe is destroyed [3]. 

We can actually considerably reduce the number of sequences that need to be checked by 
computing the rotation number of the symbol sequence using the "self-rotation number" |39j . 
The self-rotation number of symbol sequence is the same as that of the orbit, and is unchanged as 

k varies. As shown in Table 2 of [SHI the self-rotation number of ( \-)°° is ^, thus it cannot 

be one of the rotational period-four orbits — indeed it arises in a period-doubling bifurcation of 
the period-two orbit, ( — 

Proceeding in this way we identified all of the class-one rotational orbits with periods less 
than 20 40 . There are 118 such orbits; the first few of these are shown in Tabled In the 
table we chose a canonical permutation to order these sequences (we will discuss the canonical 
ordering in ^3.2)1 . 

Our numerical observations indicate that all of the class-one rotational orbits smoothly 
continue to the AI limit; we know of no proof of this observation. We observe that the codes 
for the negative and positive residue orbits differ only in one symbol, the second: indeed, this 
symbol is simply the sign of the residue, indicated by a * in the symbol sequence. 

At this point the pattern for the rotational codes that we will describe in the next section 
became clear. 



3.2 Rotational Codes 

Since the dynamics on a smooth invariant circle with irrational rotation number uj is conjugate to 
a rigid rotation with this rotation number, it is natural to construct a symbolic code for rotational 
orbits based this conjugacy. Letting 6 denote a point on the circle with unit circumference, the 
rigid rotation is 

e^9 + L0 modi et = {u;t + a}, (12) 

where a is the initial phase. Here {x} = x — [x\ denotes the fractional part of x, and [xj = 
maxpgz{p ■ p < x} is the floor function. To agree with our sign choice in Eq. we will depict 
this clS Si clockwise rotation, see Fig. 0] 

To obtain the codes for the circle map, divide the circle into two wedges 

W.iu) = i-L0,L0) 

W+iu) = [u;,l-u;], 

see Fig. 0] Note that W-{uj) is open and W+(u;) is closed, and that these wedges divide the 
circle into two parts for any < a; < ^. 

Definition 1 (Rotational (class-one) Code). A point 9 £ Ws{uj) is defined to have symbol 
s. The code for the rotational orbit 6t = ivt + a is 

st = ± when Ot G W±{oj) . (14) 

There are two distinct types of codes that we distinguish by the value of the initial phase a. 
When a = 0, the first two symbols of the code are always — h while if a S {—uj,0) then they 

are . Indeed, if a 7^ {j'^} foi' integer j, then the symbol — always appears doubled in 

the code.^ To distinguish these two cases, we call the case that a / {jw} for some j G Z a 
hyperbolic-code, and the alternative case an elliptic-code. 

■^For hyperbolic codes for rational uj this is Lemma 6.11 in |40|. 
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UJ 


Code 


k{uj) 


PD 


n 
u 

1 




-1 


3 


2 


(_+)oo 


3 


4 


1 
1 

3 


(- *+)°° 




4 


O 

4 


1 
1 

4 







0.2174036214+ 


1 
5 


(-* + + +)°^ 


/ — OV 

8 


-0.2404626622* 


2 
5 




7+5\/5 

8 


2.822983929* 


1 
6 


(- * +4^oo 


3 
4 


-0.4766507416 


1 
7 


(- * +5)'^ 


-0.8582400707 


-0.6124008240 


Z 

7 


(-* + -- + +)°° 


0.4945574340 


0.7226142786 


Q 
O 

7 


(_*_4+)oo 


2.613682637 


3.099045238 


1 
8 


(- * +6)oo 




-0.6974167690 


3 
8 




i + ^/2 


2.4246586398 


1 
9 




-0.9452647974 


-0.7546238304 


2 
9 


(-* + + -- + + +)°° 


-0.3171426657 


-0.0739574299 


4 
9 


(- * -6 + )oo 


2.762407463 


3.187557989 


1 

10 


(- * 


— i— ov 3 
8 


-0.7954212145 


3 
10 


(-* + -- + -- ++)°° 


-l+3v^ 
8 


0.8673114431 


1 

11 


(_ *+9)oo 


-0.9747995591 


-0.8258655309 


2 
11 


(-* + + + - -+^)°° 


-0.6582603930 


-0.4217831728 


3 
11 


(-* + -- + + -- + +)°° 


0.3048831897 


0.5506136016 


4 
11 


(_>,__ + __ + __ +)oo 


1.738564049 


2.225483132 


5 
11 


(_*_8 + )oo 


2.839612714 


3.215096092 


1 

12 


+10)00 


I-V3 


-0.8493921692 


5 

12 


(_*_4 + _4+)oo 


1 + ^3 


2.908515654 



t The smaller real root of lefc** - 64fc^ - 8fc^ + 1 



* The two real roots of 4096fc^ - 13312A;^ + 512fc* + 17344fc^ - 7520fc^ - 16276fc - 3251 

Table 1: Codes for class-one rotational orbits. The * represents ± for positive and negative residue orbits, 
respectively. The third column is the value Eq. at which the orbits collide with (— )°° in a rotational 

bifurcation. The column "PD" is the value of k for which the positive residue orbits undergo period-doubling. 
The elliptic i and ^ orbits actually are created earlier in saddle-center bifurcations at = 1.0 and k — 0.7063926, 
respectively. 
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Figure 4: Construction of the rotational code for oj = 7 (where 7 is the golden mean) using the wedges 
W±{u>). The orbit shown corresponds to the elliptic code with a — 0. 

Note that whenever w 7^ ^ there are at least two distmct codes. For example if we set a; = |, 

then the three possible elliptic codes are — H h and +H — , all of which represent the same 

periodic orbit with code, ( — h +)°° . Similarly all the hyperbolic codes are cyclic permutations 

of ( In the exceptional case, u = ^, there is only one period-two code, the elliptic 

code, 

It is easy to see that a hyperbolic code consists of blocks of the form |_ni-2 ^^^^ h™""^, 

where m = For example, denoting the golden mean by 7 = , then uj = 7^^ ~ 

0.381966011. Therefore 7-2 e and m = 2, so that the code is built from the blocks 

and h- For example if we choose a = —.1, we obtain 

(. .. + . + -- H + + -- + H + -...) 

For irrational uj there are infinitely many hyperbolic codes, depending upon the choice of a. By 
contrast, there is only one elliptic code, up to cyclic permutations, because any choice a = juj 
generates the same sequence of symbols. When a = the elliptic code for uj = 7"^ is 

(... + --+. - + -- H H + + H + + ...): 

in agreement with the first few points shown in Fig. |1J As previously asserted the "first" two 
symbols in this code are — +, but all other blocks have an even number of — 's. 

When UJ is rational, Def. Ogives rise to exactly two distinct codes (up to cyclic permutation). 
As we will see, these correspond to the two class-one orbits with this rotation number. 

Lemma 1. For any rational rotation number in (0, there are exactly two codes (the elliptic 
and hyperbolic codes) up to cyclic permutations. 

Proof. Set w = |, where p,q G Z are relatively prime. The periodic orbit consists of q evenly 
spaced points on the circle, and subsequent points are obtained by shifting p points around the 
circle. Partition the circle into q half-open sectors 

S, = i^-,^-^], j = 0,l,...q-l. (15) 

q q 

^This agrees with the fact that the Henon map also has only one period-2 orbit. 
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Then each Sj contains exactly one point of the orbit. Note that consists of the interior of 
the union of 2p of the Sj 

p-i 

W- = Int U Sj , 

j=-p 

(where the indices are taken mod q) and the interior of is the union of the remaining ones 

Int W+ = Int U Sj . 

j=p 

It is immediately clear that each orbit with a S Int 5*0 has the same (hyperbolic) code, since 
the points fall in the interiors of exactly the same sequence of Sj, and therefore in the same 
sequence of Wi. Similarly, for any j, the codes for orbits that have a G Int S'j are all identical. 
Finally since an orbit that starts in the interior of Sq reaches Sj after some number, t, steps, 
then the codes of orbits that start in these two sectors are the same up to cyclic permutation. 

For the elliptic case we set a = j/q for integer < j < q. Each of these orbits corresponds 
to the same set of points on the circle (the boundary points of the Sj), differing only by cyclic 
permutation. Thus their codes also differ only by cyclic permutation. 

Finally the elliptic code has only 2p — 1 points in W-{ij), while the hyperbolic code has 2p 
points. Thus these two codes differ. □ 

We choose a canonical ordering for these two codes by selecting a = for the elliptic case 
and a = — ^ for the hyperbolic case (this mimics the interleaving of the two orbits of an island 

chain). Thus the elliptic code always starts with — h and the hyperbolic with . From the 

third symbol onwards the two codes are identical. 

Lemma 2 (Properties of Rotational Codes). The canonically ordered elliptic and hyperbolic 
codes for periodic orbits differ only in their second symbol: si = — for hyperbolic and si = + 
for elliptic codes. Elliptic codes have odd parity and hyperbolic codes have even parity. 

Proof. This follows immediately from the fact that both orbits initially lie in the sector Sq-i 
defined in the previous lemma. Thus they visit exactly the same sectors in the same order. The 
only symbolic difference occurs at t = 1 where the elliptic point lies at 0i = w G W+, while the 
hyperbolic point is in the interior of Sp-i, which is in W-. The elliptic code has odd parity 
since it has exactly 2p — I points in W-. □ 

A convenient representation for rotational codes is the Farey tree. Fig. [21 The tree shows 
that the canonical hyperbolic code of a Farey daughter is simply the concatenation of the codes 
of its hyperbolic parents (putting the code from the larger rotation number on the left). This 
fact is proved in Appendix I A. 31 

As we will see in ^3.51 the wedge shaped symbol boundaries in Fig. |1J translate into a similar 
wedge for the symbol boundary in the phase space of the Henon map. However, before examining 
this, we present additional numerical evidence that lead to Def. ^ 

3.3 Numerical Observations 

While there is no proof that the codes for the rotational orbits of the Henon map are generated 
by Eq. (|14|) . we conjecture that this is the case: 

Conjecture 1. Every class-one rotational orbit of the area-preserving Henon map smoothly 
continues to an orbit of the horseshoe with code defined by Eq. il4\ )- 
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0/1 (+) 1/2 (-) 




1/5 2/7 3/8 3/7 
(..+++) (-+-++) (--+-+) (- +) 




1/6 2/9 3/11 3/10 4/11 5/13 5/12 4/9 

(--++++) (~+H h++) {--+--+H H-) ( — ] 1 H-) ( ] ] h) { 1 1 h) { ] h) ( +) 

Figure 5: Farey tree for the base 2 and i up to level three, and the corresponding codes for the hyperbolic 
rotational orbits. The code for i does not correspond to a real orbit, as there is no "hyperbolic" period-two orbit. 
The codes for the elliptic rotational orbits are obtained by flipping the second symbol. 

As noted in ^3.11 our initial observations indicate that this conjecture is true up to period 20. 
In this section we show in addition that Eq. (|14|) generates the correct codes for all rotational 
orbits up to period 99; i.e., that these orbits are born in a rotational bifurcation of elliptic fixed 
point at the value of k given by Eq. Q. To test this we must numerically continue each such orbit 
to parameter value where its residue vanishes. We observe that each hyperbolic rotational orbit 
has a negative residue near the Al-limit, and that r monotonically increases as k decreases until 
it reaches zero at a rotational or saddle-center bifurcation. Similarly elliptic rotational orbits 
have positive residues near the Al-limit, and r monotonically decreases to zero as k decreases. 
In Fig. iniwe show the scaled residue, Eq. (fTT|) . as a function k for five elliptic and hyperbolic 
class-one orbits. 




Figure 6: Left Panel: Scaled residue Eq. IIH as a function of the deviation k from k{uj) for 5 class-one rotational 
orbits. Right panel: an enlargement of the residue for the co = orbit in the vicinity of the rotational and twistless 
bifurcations. 

Thus it would appear that a quick method for detecting the rotational bifurcation would be 
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to find the parameter value, ki{uj) for the first zero of r moving away from the AI hmit. In 
Fig. [7|we plot the number of digits of agreement between ki{Lv) and k{uj) for the 1501 class-one 
rotational orbits up to period 99 as a function of rotation number. This shows that the first zero 
crossing of r corresponds to the rotational bifurcation for most of the orbits. The dashed curve 
in the right hand pane of the figure shows the cumulative distribution function for the precision 
data for the complete set of 1501 orbits. There are two groups of orbits for which the apparent 
precision is lower than 7 digits. 




to # Decimal Digits in A( to ) 

Figure 7: Left Panel: Plot of — logj^Q(fci(a;) — k{Lj)) for the 1501 class one rotational orbits with periods up to 
99. The right panel shows the distribution of digits of precision in the 1501 numerical bifurcation values (dashed 
line), and 1305 values for orbits with u! > 0.025 and either k{uj) < ^ or k{u!) > | (solid line). The latter case 
are the class-one orbits created at parameter values where the Henon map has twist and the rotation frequency 
is large enough to avoid numerical difficulties. 

One of these groups corresponds to orbits with uj G [0.29021, |], shown by the (green) 
dashed lines in Fig. [Tj This rotation number interval corresponds to the parameter interval 
^ < /c < I where the twist of the Henon map is "anomalous" . By "normal" twist we mean 
that the rotation number decreases monotonically moving away from the elliptic fixed point. 
In the anomalous interval the rotation number increases near the elliptic point before reaching 
a maximum and eventually decreasing. For orbits in this interval the first zero of the residue 
corresponds to a saddle-node bifurcation at ki{uj) < k{uj). For example, in the right panel of 
Fig. iniwe show an enlargement of the scaled residue plot for the ^ orbits. For this case, the 
first zero crossing along the curve from the AI limit occurs at /c ~ 0.7063832 for the hyperbolic 
orbit and k ~ 0.7063926 for the elliptic orbit (these points are almost indistinguishable in the 
figure). The rotational bifurcation (which is the reference value in Fig. [7j) corresponds to the 
second zero crossing along each curve; this occurs at /c(3/10) ~ 0.7135255. This explains the 
apparent loss in accuracy in the second interval (bounded by the green dashed lines) in Fig. [7| 

Our continuation method could not find the bifurcation values for orbits with w < ^ (indi- 
cated by the (red) dashed line in Fig.[7j) to more than 4 digits. This is due to the ill-conditioning 
of the equations as k ^ —1. In this region the elliptic and hyperbolic fixed points approach 
one another and the map is nearly integrable. This is reflected in the fact that periodic orbits 
very nearly lie on periodic invariant circles and so there is nearly a null direction in the matrices 
used in the continuation method. Our numerical method fails when the condition number of the 
matrices exceeds 10^^, and for low frequency rotational orbits the onset of this ill-conditioning 
occurs much further from the bifurcation value than for higher frequency orbits. This causes a 
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premature loss of precision and produces less accurate estimates of the bifurcation values. 

The solid curve in the right hand pane of Fig. [7| shows the cumulative precision data for the 
subset of orbits with lo > 0.025 and either k{u)) < ^ or k{Lo) > |. The continutation method 
produced at least 7 digits of precision for 70% of of the orbits; none of the obits had fewer than 
4 digits of precision once we excluded the orbits with anomalous twist and low frequencies. In 
every case that we investigated, the precision can be improved by carefully "hand-tuning" the 
parameters of the continuation method (in particular the parameters in the variable step-size 
routines). Thus we believe that Fig. [7| gives compelling evidence in favor of Conjecture ^ 

3.4 Symmetry 

As is well-known 41 , the area-preserving Henon map is reversible; that is, it is conjugate to its 
inverse f o S = S o f^^ ^ where 

S{x,x') = {x',x) . (16) 

Note that the reversor S is orientation reversing, det(DS) = — 1, and is an involution, S'^ = id. 
Whenever a map is reversible in this sense, it can be factored into a pair of orientation-reversing 
involutions, / = (fS) o S = R o S . The second reversor for Eq. ^ is 

R{x,x') = {x,-x' -k + x"^) . (17) 

More generally, any image of a reversor is also a reversor; the full set of involution pairs 
{{f^R, S),j G Z} form a family of reversors. Thus for example, 

T{x, x) = f^^S{x, x) = Sf{x, x) = {—x — k + x'"^, x) (18) 

is also a reversor, and since f^^R = S, we have / = ST. We will sometimes find it convenient 
to use this alternative factorization. 

An orbit O is symmetric if it is invariant under a reversor, that is if SO = O. It is easy 
to see that every symmetric periodic orbit has two points on the fixed sets of or |1J; we 
denote these, for example, as 5 = fix(S') = {(x,y) : S{x,y) = {x,y)}: 

S = {{x, x') : X = x'} , 

n = {{x, x') : 2x' = x'^ -k} . (19) 
The fixed set of T is related to Eq. (fT9|) because f^TZ = fix(/^"i?), which implies that 

f-^n = T = {{x, x') : 2x = x'^ - k} . 

The fixed sets intersect at the elliptic and hyperbolic fixed points. Following [2S1, we divide 
the lines at the elliptic fixed point, labeling the four rays with a subscript i for the ingoing half 
of the line (leading to the hyperbolic point) and subscript o for the outgoing half (leading away 
to infinity). An example is shown in Fig. IHJfor k = 2.75, where the hyperbolic fixed point is 
off-scale (to the upper right), and the elliptic fixed point is centered. Shown are a number of 
class-one invariant circles and periodic orbits that have class-two invariant circles around them. 
Each of the class-one periodic orbits is symmetric; for example, the y elliptic orbit has points 
on TZo and Si, while its hyperbolic partner has points on the other two rays. 

It was first noted by Greene that for many reversible examples the class-one, elliptic, rota- 
tional periodic orbits tend to lie on one particular symmetry ray, called the dominant ray [291 
§1.2.3]. To our knowledge this has never been proven in general, though it is known to be true 
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Figure 8: Phase space of the Henon map for k = 2.75 with a range (—1.75, 0) for x and x' . Shown are a number 
of orbits, including invariant circles around the elliptic fixed point, and five class-one island chains with rotation 
numbers f , 17 , ^ , ^ , f • Also shown are the four symmetry rays (solid curves) and their images (deished curves) 
and preimages (dash-dot curves). 

for maps similar to the standard map for large enough k [321 • -^^^ Henon map the dominant 
ray appears to be TZq for large values of k. 

However, the identity of the dominant ray depends upon the monotonicity of the twist of 
the map. For the Henon map, there is at least one anomalous domain, ^ < A; < |, where 
the rotation number is increasing -^o^ orbits in an anomalous twist region, the elliptic and 
hyperbolic orbits reverse roles, and the ray TZq contains only hyperbolic orbits. As a result there 
is no dominant symmetry line that contains all elliptic orbits when the twist is reversed. 

Under the assumption that IZo is dominant, we can easily determine which of the other 
symmetry rays contain points on the elliptic orbit, depending upon whether the numerator 
and/or the denominator of the rotation number are even or odd, see Table |5J The labeling of 
these lines is consistent with that in Fig. |H1 

The time reversal symmetry operators also have a representation for the shift map, cr, Eq. (jSJ. 
Indeed the interpretation of S is obvious if we consider its action on a configuration sequence: 
S{xt-i., xt) = {xt, xt-i). Thus as an operator on symbol sequences, 5 simply refiects them about 
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Ed 




Hd 


Hs 


odd/odd 


Ho 






So 


even/odd 




So 




Si 


odd /even 






So 


Si 



Table 2: The dominant (iJd, Hd) and subdominant(i?s, Hs) symmetry rays. Every elliptic (positive residue) 
class-one orbit (except when the twist is anomalous) is observed to have a point on Ed- Each row corresponds to 
rotation number - with p and q even or odd as indicated. 



the binary point 

S{. . . S-I.SqSi ...) = (••• SlSQ.S-l . . .) 

For symbol sequences, it is convenient to use the operator T = Sa for the second symmetry 
(instead of R), because this corresponds to a reflection about the symbol sq: 

T{. . . s_i.soSi ...) = (... si.sos^i . . .) 

Thus the elements of the fixed set of S are sequences symmetric under a refiection about the 
binary point, and those of T are sequences symmetric under a reflection about sq. Using the 
geometry of the horseshoe (see Appendix[n| and its implied ordering relation given in Lemma OOl 
we can also flnd the codes for the symmetry rays. The o rays correspond to codes whose symbols 
are "smaller" than that for (— )°^. This means that the parity of sequence sqSi . . . Sj, where sj is 
the first symbol that is not — , must be negative, i.e., there must be an odd block of — symbols 
after the binary point. Similarly the i rays correspond to codes that are "larger" than (— )°°, 
i.e. that have an even block of — 's (or none). More concretely, define k(s) to be the number of 
contiguous — symbols following the binary point:^ 

k(s) = min{j : sj = +} (20) 

then we have 



50 = {s : s-t = st~i, k{s) odd} , 

51 = {s : s-t = st-i, k{s) even} , 
To = {s : s-t = St, k{s) odd} , 

Ti = {s : s-t = St, k(s) even} . (21) 

It is easy to see that the rotational codes of ^are symmetric. For example, from Tabled the 

code s = (... !-++•••) for the hyperbolic ^ orbit has two points on symmetry lines, 

its first image cr(s) = (... + + -1- — .— -|--|--|-...)isa point on So and (7^(s) = (. . . h . H — I • • •) 

is a point on %. 

To prove this we start with a lemma about the symmetries of the rigid rotation Eq. (|12() . 
Though the reflection point for the rigid rotation is arbitrary, we choose 9 = 0, which conforms 
to the symmetry of our W± partitions, and to our selection of 9 for the canonical elliptic code. 

Lemma 3. The rigid rotation F[9) = 9 + lo mod 1 can he factored as F = ST, with reversors 
T : 9 ^ —9, and S : 9 ^ uo — 9. When uj the fixed sets of these symmetries divide into the 

^k(s) does not exist for (— )°°, as it should not, since this point is the dividing point for the four rays. 
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rays 



So ={9:9 = ^}, 5. ={9:6 = ^^}, 

% ={9:9 = 0}, % ={9 -.9 = ]^]. (22) 

Points that start on these rays generate codes with symmetries given by Eq. \21]). 

Proof. Since the W± partition is symmetric with respect to = 0, if we set = then the 
resulting code is symmetric with respect to the binary point. Since G and 9i = uj ^ W+, 
the sequence begins . — h • • •, so that k(s) is odd, agreeing with Similarly, codes generated 
by ^0 = ^ are also invariant under T, and start .+..., so they are on %. 

If we set 9q = ^, then the symmetry of the wedges (recall Fig. [IJ implies s-t = st-i, so the 
code is invariant under S. The first symbol is always sq = — ; we will show that k{s) is odd. 
Indeed, when w < |, then si = +, since 9i = ^u> < 1 — uj. Otherwise, 1 — a; < 0i < 1, so that 
si = —. But then S2 = — too, because < 92 < UJ. Indeed, a double — must occur whenever 
the orbit enters W- having skipped W+. Thus the line 9 = ^ corresponds to k{s) odd, and is 
therefore Sq. A similar argument gives the last case. □ 

Another simple lemma gives relations between the symmetry lines that imply the "dominant- 
subdominant" pairing in Table |2j 

Lemma 4. The four symmetry rays can be divided into two pairs {Ed, Eg) and {H(i,Hs), as 
given in Table \B ( depending on the parity of lo ), such that if 9 is a point on a dominant ray 

[- £— 1 -| 

then F' ~'9 is a point on its subdominant partner. Moreover the iterates of lines from different 
pairs are disjoint. 

Proof. We consider only the case to = odd/even; the others can be done similarly. First we 
prove that the stated pair of rays can be mapped into each other. That E^ is mapped into Eg is 
the statement F*[To) = % for some t. This requires that tu) = \-\- j for some integer j, which is 
equivalent to w = Thus uo is of the form odd/even, and t = f - -^o^ ^^e S pair we similarly 

have + tuj = + j, which again occurs when i = |. Since each pair maps to itself, it is 
enough to show that 7^ does not map to So, i.e. F^[To) ^ So for all t. This implies tuj ^ ^ 
mod 1. It is clear this equation does not have an integer solution for t for the odd/even case. 
Since IZ is the image of T, this verifies the pairing in Table l2j 

We find that for the odd q case T is mapped into S after 2±1 iterations. Since we use TZ in 
Table 121 this decreases the number of iterates by 1. □ 

Finally, we use these results to show that the rotational codes have the symmetries as given 
in Table 111 

Theorem 5. The codes for the rotational periodic orbits are symmetric. Moreover the canoni- 
cally ordered elliptic code is on the dominant ray To (its image is on TZo). The other symmetry 
lines occur as given in Table\^ 

Proof. The canonical code for an elliptic orbit has 9t = wt. Therefore, 9o = 0, and Lemma (HI 
implies that this point is on 7^. Since f{To) = TZo, the first image of the canonical code is on 
TZo- Thus all elliptic orbits have points on the dominant symmetry line. By Lemma^l they also 
must have a point on the subdominant line Eg. 

The canonical code for the hyperbolic | orbit has 9t = u>t — and the orbit covers the q 

evenly spaced points ,j = 0, 1, ... g — 1 on the circle. Whenever p is odd, there exists a 
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time t' such that 2j — 1 = p, and then 9^/ = £^ = ^. Thus 6t' G 5o. If (7 is odd, then there exists 
a t" such that 2j — 1 = q so that = ^ € 7^. This verifies the hne H^. The subdominant hne 
then follows from Lemma |3 □ 

In particular, this implies that the dominant symmetry line conjecture follows from the 
rotational codes conjecture. 

Corollary 6 (Dominant Symmetry Line). If Conjecture^ is true, then every elliptic class- 
one one parameter family of periodic orbits of the Henon map starting at the AI limit until 
residue zero is encountered has a point on TZq- 

The caveat in this corollary accounts for the domain of anomalous twist. The elliptic orbits 
born in this region in a rotational bifurcation have points on the H symmetry lines in Table [21 
and connect to the normal-twist hyperbolic orbits in saddle-node bifurcations when they touch 
the twistless curve. Coming from the AI limit this saddle-node bifurcation is the first encounter 
of r = 0, so the backwards running branch between k{Lo) and ki{io) is removed, because there 
the statement is not true. 

3.5 Symbol Boundary 

According to TheoremEl the boundaries of the wedges W± correspond to the image and preimage 
of a dominant symmetry line. Thus we expect that the symbol partition for class-one orbits 
will show this same structure. Indeed, this is what we observe, see Fig. 121 The figure shows 
that the wedge shaped symbol boundary leaving the elliptic fixed point (small magenta circle) 
is delineated by the image and preimage of TZq, the dominant ray. This boundary also is valid 
for the nonrotational orbits shown. 

An enlargement is shown in Fig. IIUI where there is a prominent class-two island around 
the I rotational orbit. Moving from the elliptic fixed point outward, the figure shows class-one 
periodic orbits with frequencies |g, and Notice that as expected from 

the schematic representation above, the — region is a wedge with opening angle 2a; centered on 
the dominant symmetry line. Indeed as far as we can tell the symbol boundary for all rotational 
class-one orbits is indeed formed from the image and preimage of TZq- 

Also shown in Fig. are three orbits that encircle the | elliptic orbit — orbits of class two. 
To the resolution of this figure, it appears that the wedges formed from images of TZo also provide 
the symbol boundary for these orbits. However, as we will see in the sections below, this is not 
the case. We start by constructing the rotational codes for higher class orbits. 

4 Class-Two Orbits 

As first observed by Birkhoff, a typical elliptic periodic orbit has satellite periodic orbits in its 
neigbhorhood, i.e., orbits that rotate around the elliptic orbit. Some of these new orbits will in 
turn be elliptic and therefore have satellites about them as well. We use the term class to refer 
to this hierarchy of islands-around-islands |2^. For the Henon map, the elliptic fixed point is 
defined to be class-zero, the rotational orbits of ^ are class one, and a class-two orbit is one 
that rotates around a class-one elliptic periodic orbit, see Fig. |Slfor examples. 

®It is convenient to use the Farey tree to choose periodic orbits. For k = 2.43, the rotation number of the 
eUiptic point is li; « 0.4123 < Starting with the pair of neighboring rationals | and j^, we construct the Farey 
tree for several levels to obtain our orbits. 
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Figure 9: Symbol boundary for k — 2.75 appears to coincide with the image and preimage of the dominant 
symmetry hne TZo- Shown are the same orbits as in Fig. |3 Those points with so = + are colored magenta and 
those with so = — are cyan. The solid green line is the dominant symmetry line, TZo, the dashed line is its image 
and the dot-dashed line its preimage. 



Recall that an elliptic/hyperbolic pair of class-one periodic orbits with rotation number 
a;i = 1^ is born at the elliptic fixed point when k = k{uJi), Eq. (jlUj) . After its birth, the elliptic 
class-one orbit is the center of a class-two island chain consisting of the invariant circles and 
cantori that encircle each point on the class-one orbit. The island extends to the separatrix 
formed from stable and unstable manifolds of the class-one hyperbolic orbit, as schematically 
shown in Fig. 1111 The class-two orbits have rotation numbers that are measured relative to 
fii . When the linearized rotation number of the class-one orbit — as determined by its residue, 
Eq. (O or Eq. Q — passes through the rational frequency P2/Q2, an elliptic/hyperbolic pair of 
class-two orbits is created in a g2-tupling bifurcation of /'^. We denote the rotation number of 
these orbits by wi : 102] they have period Q = qiq2- Three such class-two orbits are shown in 
Fig. El with rotation numbers | : | : ^ and | : ^• 

For the case of normal twist, the direction of rotation typically reverses with each increment 
of class. This occurs because in an any island the rotation number is typically a maximum 
at the elliptic orbit forming the island center and monotonically decreases to zero approaching 
the separatrix. In particular this implies that the inner separatrix of an island advances with 
respect to the center of the island and the outer one retreats, see Fig. ^2 Continuity then 
implies that if the class-one invariant circles have clockwise dynamics, the class-two circles have 
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-1.2 -0.8 X -0-4 

Figure 10: Class-one symbol boundary for k « 2.43 (e = 0.54) coincides with the image and preimage of the 
dominant symmetry line IZo- Shown are 13 class-one rotational orbits about the elliptic fixed point, as well 
as three class- two orbits encircling the | elliptic orbit. Elliptic (positive residue) orbits are shown as dots and 
hyperbolic as crosses. 

counterclockwise dynamics. 

Our numerical observations imply that the code for a class-two — : — orbit is constructed 

^ J qi q2 

essentially by repeating the class-one ^ code q2 times, because most of the class-two points are 
deeply buried inside the wedges shown in Fig. EJthat determine the class-one symbols. However, 
if the class-one sequence were merely repeated q2 times, the orbit would not have the correct 
minimal period; consequently some of the symbols must be flipped to obtain the correct code. 
Indeed, there are two islands that intersect the class-one symbol boundary (recall the | island 
in Fig. nr^ J Such islands are termed ambiguous because they straddle the symbol boundary; 
the determination of the symbols of points in these special islands requires more analysis. 

4.1 Class-Two Codes 

Schematically, a class-one orbit with rotation number loi = ^ can be represented as qi equally 
spaced points on a circle with the dynamics Eq. (|12() . In a similar vein, a class-two orbit with 
rotation number 2i : £| schematically corresponds to points on a set of circles enclosing each 
point of the class one orbit. To compute these points choose a pair of radii, ri, r2 to represent 
the sizes of the islands and define 

x + iy = e"2^*^i (ri + rse^"^^^ . (23) 
^Except for the case oJi = i, where there is only one. 
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Figure 11: Schematic class hierarchy: class-one invariant circles (cyan) and a class-two j island chain (red). 
Since the rotation number of the class-one circles decreases with radius, the image of the horizontal (green) line 
shears as shown. Continuity implies the direction of rotation reverses as the class is incremented. Also sketched 
are the elliptic (green) and hyperbolic (brown) symmetry rays for class-one. 

Here r2 chosen sufficiently smaller than ri so that none of the class- two islands intersect. This 
gives an epicycle view of the orbit. The angles (6*1, 62), representing the island structure, evolve 
according to the two rotation numbers of the orbit. The schematic dynamics is expressed by 
using a mixed basis for time: t = ti + qit2 where ti = t mod qi so that ti G [0,qi); this is 
written in clock-like notation as® 



Iteration of this pair of maps produces a set of Q = qiq2 points on the epicycle. 

The signs in the exponentials in Eq. ()23() are reversed to account for the reversal of rotation 
direction that occurs when moving between the classes. 

Recall that the canonical elliptic class-one code is obtained by setting ai = and the 
hyperbolic code is obtained when ai = — l/2gi. Since class- two orbits rotate around elliptic 
class-one orbits, we set ai = 0; the phase 02 differentiates between the elliptic and hyperbolic 
class-two orbits. The canonical elliptic class-two code has 02 = 0, and the canonical hyperbolic 
code has 02 = 

The epicycles for a; = | : | are shown in Fig. 1121 The points of the elliptic class-one orbit 
are represented by filled circles evenly spaced along the large circle. Surrounding each of the qi 
points of this orbit is a smaller circle containing q2 points of the class-two periodic orbit. By 
convention, we measure 62 in a corotating coordinate system, so that points with 62 = are 

^Mixed base systems like the above are well known for time measurements, e.g. setting lo = : : : 
^ : i : ^ corresponding to tenth of seconds, seconds, minutes, hours, weekdays, weeks, so that in our notation 
t = : 1 : 30 : 19 : 4 : 5 is the 4th day in the 5th week, at 19:30 hours plus 1.0 seconds. Instead of the American 
system uses : |, where the 2 is represented as AM or PM. Our lives are simpler because all numerators are 
1 for an actual clock, hence dj is not needed; they would have been even more so if the French revolution had 
succeeded in promulgating a decimal time system. 



t = ti:t2. 



With this definition, the angles evolve as 




(24) 
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located on radial rays from the origin. Notice that a gi-fold iteration maps each small circle to 
itself, shifting its points counterclockwise by uJ2- 




Figure 12: Epicycle view of a class-two orbit for lj = 2 : |. The picture is generated by Eq. 11231 with ri = 1 
and r2 = 0.1. The angle 6i is measured clockwise, and 62 counterclockwise. The elliptic class-one | orbit is 
denoted with solid circles, and the class-two orbits are denoted by open circles (elliptic), and crosses (hyperbolic). 
Points with symbol — are colored blue and those with + are colored red. The integers next to each elliptic point 
denote the time t along the orbit. 

To generate the symbol sequence for a class-two orbit, the epicycle picture must be parti- 
tioned into regions corresponding to the + and — symbols. Our observations indicate that this 
can be done by constructing a pair of "wedges" W±{ijj) whose boundaries are piecewise linear. 
The boundaries start as rays from the origin at angles Oi = ^uji — these same rays formed the 
class-one boundary. The rays end in the center of the two ambiguous class-two islands corre- 
sponding to ti = ±1; it is not immediately clear how to continue the symbol boundary in these 
islands. A cursory look at Fig. might indicate that the symbol boundary continues straight 
through the ambiguous islands. However, our numerical observations show that the two rays 
actually bend at the centers of the = ±1 ambiguous islands by angles 62 = ±1^2, respectively. 
This bending is too small to see in Fig. Ilfll 

The regions : UJ2) bounded by these rays correspond to the it symbols. Just as for the 

class-one case, the region containing 61 = Q corresponds to St = — . The only remaining point to 
understand is the symbol for points on the boundary. In contrast to class-one, our observations 
show that W- is closed and Wj^ is open. 

Definition 2 (Class-Two Code), t = ti : t2 determines a point z in the plane by \2^) and 
\24\) - If z £ Ws{lji : L02) then the tth symbol is s. The region Ws{uJi '■ UJ2) is bounded by a 
continuous piecewise linear system of segments. At the origin two segments of length ri emerge 
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at angles corresponding to ti = ±1. Their endpoints are connected to segments of length r2 with 
angles corresponding to ti : t2 = ±1 : ±1, see, e.g. Fia. liM The region W-.{uJi : UJ2) is the one 
that contains t = : and is closed, while W^{uJi ■ UJ2) is its complement. 

It is not difficult to translate the epicycle picture into a simple algorithm that generates the 
codes for any io = — : —. 

Continuing the example of Fig. 1121 the orbits are coded so that the points in the red region 
have symbol + and those in the blue region have symbol — . The four points that fall on the 
symbol boundary (for example t = 8 = 1 : 1 and t = 62 = 6 : 1 for the elliptic orbit) have 
the symbol — . Note that points in the non-ambiguous islands, corresponding to ti = 0,2,3,4,5 
have the same code as the elliptic class-one y orbit. Thus many of the symbols of the class-two 
orbit simply repeat the class-one code. 

To see this more clearly, denote the class one code by the row vector si. Since most of 
the points on the class-two orbit are determined by this code, we first replicate this code Qc 
times; however, instead of viewing this as a vector of length qiq2 we reshape the code as a table 
with si as its rows. The columns of this table are labeled by ti, and the rows by t2. In this 
way each ambiguous island corresponds to a column of the table, i.e. the columns ti = 1 and 
ti = — 1 = gi — 1. The codes in the ambiguous columns must now be replaced by the appropriate 
code for L02 using the class-two wedges. 

As an example consider again the symbol sequence the I : I orbit. The Table El hsts the 
symbol sequence where the 63 symbols are arranged in 9 rows of 7 entries each. The ambiguous 
islands correspond to the columns ti = 1 and ti = 6 = — 1 mod 7. The 92 (^Zi — 2) symbols 
in the remaining columns are determined by g2-fold repetition of the elliptic class-one code for 
a;i = 2/7, {- + + -- + +r . 
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Table 3: Elliptic class-two code for the | : | orbit. The ambiguous columns are ti = 1 
and ti = —1. The hyperbolic code is the same with the exception of the boxed symbols, at 
i = 8 = 1 : 1 and t = 50 = 1 : 7, which are ffipped. The right table compares the code for the 
elliptic I class-one orbit with the two ambiguous columns in the | : | orbit . 

The symbols in the ambiguous columns ti = ±1 are generated with respect to a wedge 
formed by the boundary in the class-two island. The codes for the hyperbolic orbit are obtained 
by a shift of 02 = "252' ^ clockwise rotation) as shown in Fig. 1121 this corresponds to 

flipping the two boxed symbols in Table |31 As we will see below, when there are two ambiguous 
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islands, exactly two symbols flip relative to the elliptic case. 

Epicyclic pictures for additional class-two orbits are shown in Fig. 1131 The codes are easily 
obtained from the figures by reading off the symbols according to the time t, using the colors 
shown. For example, the code for the elliptic ^ : y orbit is the 7 fold repetition of the ^ elliptic 
code, — h^, modifying only two symbols in the columns = ±1, by setting sio = sq2 = — • The 
hyperbolic code is obtained from the elliptic one by flipping sio to +, and S55 to — . 

As LOi approaches |, the class-one + wedge shrinks to a single point 9i = ^. Orbits with 
uji = ^ are special because the two boundaries of the class-one wedge coincide so that there is 
only one ambiguous island with ti = +1 = — 1. Thus for example, the code for the ^ : | elliptic 
orbit shown in the last panel of Fig. ^Jis ( — I — i4-joo^ while its hyperbolic partner has code 
(- + -^^+)°°. In this case the ambiguous column of the elliptic orbit contains H — ^ while the 
unambiguous column simply contains — ^. 

If UJ2 approaches ^, as in the first two panels of Fig. ^1 the open class- two wedges grow 
encompassing all but ^2 = ^- The extreme case of this, uj2 = ^ has only one s = — point in the 
ambiguous islands; for example, the code for the elliptic | : ^ orbit is ( — \-^—)°°. 

The well-known "doubling" substitution rule 

> - + and + ^ (25) 

(see Appendix El can be used to generate the codes for class- two orbits of the form ^ : u) from 
those of the class-one orbit with frequency uj. For example, applying this substitution rule to 
the I code, — gives the same code as in Fig. 1131 That this is true in general follows because 
the class-two wedge for s = + is 82 G {—uj,uj), which is the same as the class-one wedge for 
s = —. Hence each symbol flips to create the ambiguous column; the first, unambiguous column 
is all — . 



4.2 Numerical Observations 

In this section we present numerical observations that support the class-two code construction 
for the area-preserving Henon map. As for the class-one case in ^3.31 we checked that the 
class-two orbits defined by their codes at the AI limit collide with the appropriate class-one 
orbit at the parameter value where the class-one orbit has the expected residue. We restricted 
our attention to class-two orbits with qi and q2 < 13. Fig. shows the number of digits of 
precision in the bifurcation value for the class-two rotational orbits with normal twist. ^ When 
= 1/2 there is only a single ambiguous class-two island. In order to verify the class-two code 
for this case we compared the numerical bifurcation parameters for 236 class-two orbits of the 
form i : uj2 where §2 < 39 to the exact value jSO] 

: 002) = 3 + sin(7rix'2)^ . 

Outside of the interval below a;2 = 1/3 where the twist vanishes at class- two the computed values 
of the bifurcation parameter were accurate (to at least 8 significant digits) when compared to 
the exact value. 

The results we have obtained give us confidence to formulate the following 

Conjecture 2. The AI code for the class-two orbits of the Henon map is defined by the wedges 
W±{uJi : UJ2) given in j[ 

^There are only eight orbits with fewer than four digits of precision: ^ : : : • • "^'"^ • 

Ta ■ T3' ■ These are probably related to anomalous twist either near i or near j. 
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# Decimal Digits in k(a^ : co^) 



Figure 14: Number of correct digits in parameter value for 727 class-two orbits with qi < 13 and q2 < 13. We 
excluded frequencies P2/?2 = 1/3, 3/10 where the class-one twist is anomalous to avoid numerical difflculties with 
multiple bifurcations. 



4.3 Parity of Class-Two Codes 

We now show that the elhptic and hyperbohc class-two codes differ by two symbols, but they 
have the same parity. 

Lemma 7 (Properties of class-2 codes). Assume uj ^ : ^. Then the elliptic and hyperbolic 
class-2 codes differ by flipping two symbols, the symbol at t = 1 : 1 and that at 

i = -l:| 

t = l:'-^ 

where j is the smallest non-negative integer for which an integer solution for t2 is obtained. The 
parity of both orbits is equal to the parity of q2 ■ 

Proof. We start by considering the ambiguous island, ti = +1. By definition, the hyperbolic 
and elliptic codes differ only when 102^2 and Lx)2t2 — ^ are in different wedges. Let Sj be the 
sectors defined in Eq. H15|) relative to q2 for 62. When q2 is even the situation is similar to 
class-one, because the the wedge boundaries in the ambiguous island coincide with boundaries 
of some sectors Sj. Accordingly only the symbol t = 1 : 1 changes. Alternatively when q2 is 
odd, the other wedge boundary at ^ is not part of the elliptic orbit since u>2t2 ^ ^ for all ^2- 
Instead the hyperbolic orbit hits this boundary when (^2^2 — ^ = ^ + J ? flipping the symbol at 
this location. This gives the above formula. In order to show that this is the only other change 
with t = 1 : t2 the argument of Lemma ^ can be repeated with sectors Sj of half the size. As a 
result for all 02 € S-i the elliptic code is obtained, while for 02 G -S'_2 the hyperbolic code is 
obtained. 

Now we consider the ambiguous island with t = —1 : t2- The wedge in this island is a 
reflection of that in the previous ambiguous island. However, since the rotation direction is not 
reversed, the sectors Sj (or Sj) are not reflected. This implies that whenever an open sector end 
coincides with a closed wedge end for t = -|-1 : t2; after reflection of the wedge for t = —1 : t2 
both endpoints will be open, and vice versa. Since a symbol flip only occurs when different types 



q2 even 
q2 odd 
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of endpoints coincide only one of the symbols at t2 = ±1 flips, never both. In particular when 
q2 is odd both symbol flips occur for t2 = +1, while for even q2 the symbol half way around the 
orbit with 002^2 = f + J Qips for t2 = — 1- This gives the first part of the lemma. 

When q2 is even the parity of the unambiguous symbols (ti 7^ ±1) is even, because their 
number per island is multiplied by q2- The parity of the ambiguous symbols is even in the elliptic 
case, because the symbol sequences are cyclic permutations of each other: Starting with 02 = 
at = 1 : 1 gives the same symbol sequence as starting with a2 = ^ — UJ2 at t = —1 : ^ (apply 
a rotation!). But 02 is a point on the orbit for even q2, hence the elliptic ambiguous columns 
have the same numbers of signs. Therefore their parity is even, and the whole sequence has even 
parity when q2 is even. 

When q2 is odd the parity of the unambiguous symbols is odd, because the number per island 
is odd and they are repeated an odd number of times. The number per island is odd because 
the elliptic class-one orbits have odd parity, and the two deleted symbols in the ambiguous 
column are + signs in the elliptic case. The number of signs in the ambiguous columns is the 
same because the two are time reversals of each other: upon reflection of the wedge and the 
orientation the same sequences are generated. So their parity is even, and the total parity is 
odd. □ 

Recall from Lemma |^ that the elliptic/hyperbolic class-one orbits have opposite parities and 
this corresponds to their opposite residue signs. One important consequence of the previous 
lemma is that at class two the elliptic/hyperbolic partners have the same parity and therefore, 
according to Eq. Q, will have residues with the same sign at the AI limit. Nevertheless they 
are initially born with opposite residues, as this is the normal form of the g2-tupling bifurcation 
||-i4j . According to Eq. Q even parity implies a negative residue at the AI limit. Therefore, 
when q2 is even the residue of the (initially) elliptic orbit, which is (initially) positive, must cross 
r = at least once as k increases from the bifurcation value in the approach to the AI limit. 
Similarly, when §2 is odd the residue of the (initially) hyperbolic orbit must cross r = for some 
k larger than the bifurcation value. 

We call bifurcations related to these additional zeros in the residue function secondary bi- 
furcations. In particular when r = for k larger than the parameter value of the rotational 
bifurcation, a secondary pitchfork bifurcation occurs (super- or subcritical, depending on the 
parity of the code), and when r = 1 a secondary period-doubling bifurcation occurs (again super- 
or subcritical, depending on the parity). Between these two parameter values, a complete se- 
quence of secondary rotational bifurcations must take place. An example of these bifurcations 
is shown in Fig. El 

We will not attempt to completely describe the symbol sequences associated with these 
secondary bifurcations, but only offer an empirical rule based on our numerical observations: 

Conjecture 3 (Pitchfork class- two). The pitchfork children of class-two orbits only differ in 
the four symbols at the intersection of the ambiguous columns and rows with ti and t2 given in 

Lemma^ When q2 is even these four symbols ordered by t are — I for the (initially) elliptic 

orbit and + H h for the (initially) hyperbolic orbit. The secondary bifurcation of the elliptic 

orbit is 

- + .pf(- + -+,- + +-). 

When q2 is odd they are — h H — for the elliptic and + H for the hyperbolic orbit. The 

secondary bifurcation of the hyperbolic orbit is 

+ + > pf(+ + +-, + + -+) . 
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Pitchfork child 



10° lo' , 10' lo' lO'' 10 ' lo" lo' , 10' lo' 10* 

k k 

Figure 15: Scaled residue p, Eq. Ijll^ . for typical class-two rotational orbits and their pitchfork children with 
q2 odd (left panel) and even (right panel). The (initially) stable class-two orbit is rendered in blue, the (initially) 
unstable orbit in red, and one of pitchfork children in green. The dashed curve is the period-doubling curve r = 1. 
The secondary pitchfork bifurcation occurs at the rightmost zero crossing on the (initially) stable orbit when 52 
is even and the initially unstable orbits when §2 is odd. 

It is not hard to show that the four ambiguous symbols in the parent sequence must be 

=F H as given in Conjecture 13 these are the symbols corresponding to the points on the 

appropriate class- two wedge boundaries. 

The rule given in Conjecture El generates the symbol sequences of the daughter orbits created 
in a secondary pitchfork bifurcation. These are the orbits whose residues are rendered in green 
in Fig. 1151 We examined 105 class-two orbits in the Henon map with frequencies of the form 
J : 2 where j = 3, 4, 5 and 3 < g < 15. The parameter value {k) where the pitchfork children 

predicted by ConjectureElwere born agreed (to at least eight significant figures^'') with the value 
where the residue of the class-two parent orbit crossed zero in the "secondary" bifurcation. These 
results suggest that the conjecture is probably valid for the Henon map. 

4.4 Symmetries 

Numerical observations suggest that the concept of "dominant" symmetry line can be extended 
to class-2 orbits (and in fact to class-c with c > 2, see below); for example, for each wi, there 
appears to be a particular symmetry ray that contains all of the elliptic class-two orbits loi : L02 
pS] . Recall from Table HI that there are two rays, and Eg-, whose identity depends upon the 
parity of the numerator and denominator in wi, that contain elliptic class-one orbits. Since each 
class- two orbit rotates around the points on a class-one orbit, they should have points on the 
two E rays. These two rays are divided at the class-one elliptic points, giving four new rays. 
As before, we use i and o to denote "inward" and "outward" halves of the rays; the subscript 
i corresponds to the ray that starts at the elliptic class-c orbit and heads in toward the elliptic 
class-(c— 1) orbit, while o is the ray that heads out from the lower class orbit. Thus, for example, 

(c) (c) (c) 

Es is divided into Eso and Esi- We denote the four symmetry rays at class-c by E"-^' , Ef , H'^' , 

^"With the exception of a six orbits: | : ■ yi, j '■ j^, j ■ j^y j ■ Ja' ^ ■ which numerical difficulties 

gave less than six digits of precision. 
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(c) 

and Hs ■ The observation in |25] is that 



-'^w — ^so > 



(26) 



that is, the outward half of the class-(c — 1) subdominant ray becomes the dominant ray for 

2 . P2 
5 ' <?2 



class-c, see Fig. 1161 This is consistent with Fig. Iin( as all of the elliptic, class-two 1 : ^ orbits 



have a point on the outward half of E^^^ = So- 






Figure 16: Sketch of the dominant symmetry line at class-two. Shown are points on a period 5 elliptic class-one 
orbit (cyan dots), and several class- two orbits (red on the symmetric islands) All of the elliptic class- two orbits 



line up on the outward half of the subdominant class-one ray Elo 



(1) 



In this section we verify the dominant symmetry conjecture for the class-two codes. We also 

(2) 

obtain a recipe to describe the boundary of the class-two wedge in terms of iterates of the , 
so that the wedge can be computed in the Henon map and compared to the prediction obtained 
from numerical continuation. 

That the class-two elliptic codes are symmetric is a simple consequence of the symmetry of 
the wedges W± in the two epicycle picture. 

Lemma 8 (Symmetries of class-two codes). Let s be the code for the elliptic class-two orbit 
with rotation number uj = ^ : Then s has an image on the dominant symmetry line E^ = 



E^^^ . In addition, when q2 is odd, s has an image on the subdominant line E^"^ 



E 



(1) 



no. 



and when q2 is even on E, 



(1) 



Here are the lines in Table\^ 



Proof. Let s denote the canonical elliptic symbol sequence, i.e. starting from the point 6 = (0, 0). 
Then if = 1 : 1 = 1 -|- qi, the point (T*+s is on the symbol boundary since it corresponds to 
9 = {uji,uj2)- Similarly setting t_ = — 1 : — 1 = {qi — 1) : {q2 — 1) = qiq2 — 1 = —1, then cr*-s 
is also on the symbol boundary, and corresponds to 6* = {—uJi, —i^2)- Since the class-two wedge 
W- is closed, these two points have the same first symbols so = — , and by symmetry of the 
wedges in i— > —9, the two sequences must be the same when reversed; consequently. 



Ta^+s 



a 
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Repeatedly composing this with a, and noting that aT = Ta , we find that when + t_ is 
even, a°'s £ T where a = = Alternatively when + 1_ is odd, since aT = S we find 

a'^s e S where b = t++L^ = 21+1. 

To decide which half of the symmetry line these points are on, recall from Eq. H20|) that this is 
determined by k{s), the number of consecutive — symbols starting at sq. Consider first qi even, 
where the symmetric point is 0"°s. Since t = a = ^ : 0, then 6 = (^^^,0) = (^,0) = (^,0) 
because pi is necessarily odd. Thus the first symbol of this sequence is sq = +, and so by 
Eq. (|22j) it is on %. When qi is odd, the symmetric point is cj^s, with t = : 0, and so 

9 = + ^,0). When pi is even, this is on the line 6i = with sq = — ; following the 
argument in Lemma this implies the point belongs to So- Similarly when pi is odd, the point 
is in Si. 

Summarizing we obtain the table: 





ti:t2 


01 


O2 




odd/odd 


91 + 1. Q 

2 


1 _l_ a;i 

2 2 





Si 


even/odd 


91+1. Q 


2 





So 


odd/even 


21-0 


1 
2 





% 



Since this symmetry is independent of 0J2, we are justified in calling it a "dominant" symmetry 
line. 

The subdominant symmetry line can by found by noting that since s has period Q = qiq2, 
we also have 

rcr*+S = cr*-+'3s 

Iterating this relation we find that if qi{q2 + 1) is even, then a'^s S T when c = = 
3li3^. Otherwise, a'^s E S where d = M£2+i)±l. 

The choice of ray for the subdominant symmetry depends upon the parity of both uji and 
UJ2- When q2 is odd, then the second line is at a'^s. Indeed t = c = : and so = 

independently of uji, and ^2 = ^ + For any p2 we have a'^s = . — h • • ., so that this point 
is on To. 

We now consider 52 even. If qi is also even then the symmetry still occurs at t = c = ^ : 
which gives 6 = (^,^) = (|,^) since both p's are odd. If < |, then So = +, and the 
symmetry is 7^. A similar argument gives the same result for the frequency ^ case. The final 
case has qi odd and q2 even, for which the symmetry corresponds to t = d = ; gg that 
9 = + ^). When pi is even, 9i = so sq = —. This point is on So, as follows from an 
argument similar to that in Lemma |21 Finally if pi is even k{s) is even, so the symmetry is Si. 

Summarizing, we obtain a table for the subdominant line: 





UJ2 


ti:t2 


9i 


92 




any 


even/odd 


n.92+1 
u. 2 





L02 
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To 


any 


odd/odd 


Q.92 + 1 

u. 2 





1 I 1^2 

2 2 


To 


odd/odd 


odd/even 


gi+1.92 
2 • 2 


2 2 


1 
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S^ 


even/odd 


odd/even 


gi+i . 92 
2 • 2 
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1 
2 


So 


odd/even 


odd/even 


91 .92 
2 • 2 


1 
2 


1 
2 


% 



□ 

So far we have shown that the class- two codes are consistent with Eq. (|26() . but we have not 
yet shown how to divide E'^^ into halves. 
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Theorem 9 (Class 2 Dominant Symmetry Ray). The dominant symmetry ray for class- 
two codes is the outward half of the subdominant ray at class-one, e''^^ = si]) , thus verifying 
Eq. i2b}) for class-two. 

Proof. The division of the symmetry hnes for the shift a into rays at the class-one orbit uses 
the ordering relation discussed in Appendix [O As shown in Fig. 1231 the outer halves of 5o and 
To correspond to codes that are "less" than (under this special ordering) those of the class-one 
orbit, while the outer halves of Si and % correspond to codes that are "greater" than those of 
the class-one orbit. The canonical elliptic class- two code is the same as that of the class-one 
orbit up to time t = 1 : 1 = 1 -\- qi, where it first hits the class- two wedge boundary. Moreover, 
according to Lemma|Sland LemmalU if s is the canonical elliptic code, then its [^^g^Jth iterate 

17(1) 

IS on Eg 

Thus the relative ordering of the class-one and class-two codes, following Lemma OOl depends 
upon the parity of the symbols of the class-one orbit starting halfway around the orbit, and going 
up to the first differing symbol Sg^+i = si, i.e. on the sign of 

p = vr(S|^9i+ij . . . Sgi-iSQSi) , 

In particular, we must show that p = (— 1)^^. This follows because if pi is even, Ei^^ = So, and 
so the code for the class-one orbit should be greater than that of the class-two orbit, while if pi 
is odd, then the symmetry lines are either Si and TZi, so the class-one code should be less than 
that for class-two. 

As noted in Lemma |21 the class-one elliptic code has 2pi — 1 symbols that are — . Moreover, 
the canonical elliptic code has symmetry T, which means it is unchanged under refiection about 
So = — • This implies that the first half of the orbit siS2 ■ • • s. g^-i . is the same as the second 

L 2 J 

half s 1 91+1 1 . . . Sq^-i written backward. Consequently they each have pi — 1 symbols that are 

L 2 J 

— . Finally, the sequence S:qi+i , . . . Sg^-iSoSi has one additional minus sign, making pi in total; 

L 2 J 

therefore, p = (— l)*'^. □ 



4.5 Symbol Boundary 

In order to connect the above findings to the Henon map we now describe the wedge boundaries 
in terms of iterates of symmetry lines. Inside the class-one island the wedge boundary is simply 
given by the forward and backward iterate of the dominant symmetry ray TZo, as discussed in 
H3.5I From the schematic picture it would seem as if the higher class wedges are also constructed 
from iterates of TZo- This would, however, contradict the observation that the dominant symme- 
try line is different at the next class l25|. Moreover, we expect that a symbol boundary should 
be associated with a line containing all the elliptic points. Thus, as we saw in Lemma|Hlit should 
be associated with the subdominant elliptic line. 

Therefore we now show that appropriate iterates of TZo in the schematic picture can in fact 
be identified with certain iterates of the subdominant lines as given by Eq. (|26j) . 

Lemma 10 (Boundary Symmetry Lines). The boundaries of the class-two wedge in the 
ambiguous islands ti = ±1 for the class-two orbit iv = pi/qi : P2/Q2 (',re given by 

• the [^^^J + ti iterate of the subdominant symmetry ray E^^^ and 

• the —[^^^1 iterate of symmetry ray e'^^J . 
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Proof. Recall the proof of Lemma 0J The schematic boundary of the class-two wedge has two 
parts. The inner part consists of the class-one rays F^^{TZo) with ti = ±1. The first part of the 
lemma is the statement that this ray is an image of E^^^ = So- For example, when toi = e/o, 
then we require that there exists an integer n such that 

^1 = ^- ^*n^o)=i^"(5o) ^i(ti + l)+j = ^+^in 

odd 2 

for some integer j. Solving for uji gives 2j'/(2n — 2ti — 1), which is of the form even/ odd. When 
ti = 1 the explicit solution is n = (gi -|- 3)/2, and when ti = —1 it is n = (gi — l)/2. Note that 
replacing Sq by 5j or TZi leads to a contradiction. Hence So is the unique symmetry ray that 
can be mapped to TZo for uji = even/ odd. In a similar way the other cases lead to 

^1 = ^: F'^{no)=F'\S,) a;i(ii + l)+j = i±^+a;in 

which gives the same values of n as before. Finally, 

u, = ^: F'-{no)=F^{ni) a;i(ti + l)+j = i+o.i(n+l) 

even 2 

which gives n = ^ + ^i- Combining theses expressions gives the first form. 

The second ray of the class-two wedge is formed by taking the o part of the ray just con- 
structed (i.e. the part that goes beyond radius ri), and turning it by ±cj2 for t = ±1, respectively. 
Turning by iba;2 is achieved by iterating one complete cycle of class-one, i.e. iterating itgi times. 
For ti = — 1 the number of iterates is thus {qi — l)/2 — qi for odd qi and qi/2 — 1 — qi for even 
qi. Combining both cases gives — \{qi + l)/2], and the result follows. □ 

As a result of the previous Lemma the number of iterates needed to construct the class-two 
wedge only depends on qi, the period of the class-one orbit. 

It seems pointless to construct the class-two wedge from iterates of the new dominant sym- 
metry line instead of from iterates of TZo, since the two agree in the schematic picture. It turns 
out, however, that the former gives the correct wedges for the Henon map. In the schematic 
pictures the two constructions coincide because the map has no twist. At the birth of the class- 
two island (when the rotation number of the class-one orbit is 1^25 and the rotation number of 
the class-two orbit is 0) the schematic picture is an accurate description of the wedge-geometry. 
Away from this bifurcation the symbol boundary differs from that in the schematic picture, and 
the wedge needs to be constructed from the dominant lines of class 2. 

In Fig. I17( we show the symbol boundary for k = 2.684, where the | island is predominant 
(the same phase space is shown in Fig. I18p . Here the class-one orbits are shown with colors cyan 
and magenta for s = it, respectively. The boundary between these colors is given by the image 
and preimage of TZo as before. Class-two orbits with a; = | : ^ are shown with colors blue 
and red for s = it. Note that the class-one partition boundary does not work for these orbits. 
According to Lemma |1] the boundaries in the ti = —1 ambiguous island are given by segments 
of f^{So) and /^^(5o); these are shown as the dashed (brown) curves, and clearly delineate the 
symbol boundary in Fig. El for the class-two orbits. Similarly segments of f^{So) and f^^{So) 
form the class-two boundary in the ti = 1 ambiguous island. 

There is apparently a jump in the symbol boundary across the separatrix of the | orbit. We 
are unable to resolve the behavior near the separatrix numerically as the period of the orbits 
near the separatrix becomes too large to accurately find the bifurcation values numerically. 

The ability to describe the wedge boundary that is obtained by continuation from the AI 
limit in terms of iterates of the symmetry lead to Conjecture |2I 
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Figure 17: Class-one (cyan and magenta) and class-two (blue and red) rotational orbits for k = 2.684. The 
— symbols correspond to the cyan and blue points, and the -I- to magenta and red. Also shown are the t = ±1 
iterates of the symmetry line TZo (green) that form the class-one boundary, and the iterates t = —1, 4 and t = —3, 2 
of So (brown) that form the class-two boundary. 



5 Class-c Codes 

Orbits of arbitrarily high classes can also be found in area-preserving maps. For example, Fig. 1181 
shows successive enlargements of the phase space of the Henon map and exhibits orbits up to 
class four. The enlargements in this figure are not done with the same multiplier; indeed, though 
there are special parameter values for which the class hierarchy exhibits self-similarity ^25^, in 
general it does not. While there are typically islands of every class, their shape and structure 
changes with class. 

A class-c orbit rotates around a class- (c—1) elliptic periodic orbit with some definite rotation 
number. We denote the rotation number by a string 

u; = ui : a;2 : . . . : . (27) 

Each of the numbers uJi = ^ is necessarily rational, except possibly uJc- If the class-c orbit is 
periodic, it has period Q = 11^=1 The rotation number tOc+i about a class-c periodic orbit 
is obtained by iterating and considering rotations about one point on the class-c orbit — if 
an orbit returns to itself after qc+i iterations of and has undergone Pc+i full rotations, then 

Generalizing the epicycle picture in Fig. 1121 we can construct the symbol sequences for 
periodic orbits of arbitrary class. As before, time is written in the mixed Qi basis: 

t = ti + qi{t2 + q2{t3 + ... + Qc-itc)) = h : t2 : . . . : ] (28) 
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Figure 18: Phase space of the Henon map for k = 2.684. Enlargements show higher class orbits, ultimately 
focusing on the | : ^ : ^ : orbit. The class-one invariant circles (cyan) encircle the elliptic fixed point, class- 
two (red) the | elliptic orbit (red). Shown in the enlargements are class-three circles (green), class-four (orange), 
and finally a class- five chain (mauve). The elliptic class- two | : ^ orbit has a point on the subdominant elliptic 
line So for the elliptic | orbit, as implied by dominant symmetry line conjecture. 



where the digits ti are taken modulo qi. Associated with each t there is a point 9{t) G T'^ defined 
by the mapping 

t 6{t) = {uJiti,UJ2t2, ■■■ , UJctc + Qc) • (29) 

Here we have set aj = for j = 1,2, ... c — 1; this represents the selection of elliptic orbits up 
through class-(c — 1). The phase Oc selects an elliptic or hyperbolic class-c orbit. A class-j island 
consists of the set of points 0{t) with a fixed choice of ti : t2 ■ . . . ■ tj-i. Incrementing tj moves 
the orbit around the island with frequency ooj, while incrementing 9j by ^ moves sequentially 
through the points in the island. 

The epicycle picture is obtained by generalizing Eq. (jSHJ- Define a set of radii, Vi, i = 1, . . . c, 
that decrease rapidly enough so that the set of circles defined by 

a; + ,y = (n + e^-^^ (r2 + . . . r,e(-i)^2..9.^^ ^ (30) 

do not overlap. As before, the direction of rotation reverses with each class. An example epicycle 
picture for class three is shown in Fig. 1191 

The definition of the class-c code is complete once we specify the boundary between the 
"wedges" W± that define the symbols = it. As for class- two, our numerical observations 
indicate that the boundary consists of another wedge. The rays that bound the wedge leave the 
origin with angles 9i = ztoii . These extend to the center of the two class- two ambiguous islands 
at ti = ±1. The rays turn by angles 02 = ±1^25 and continue to the center of two class-three 
ambiguous islands at t = 1 : 1 and t = —1 : —1, respectively. The rays again turn by angles 
63 = ±UJ3, continuing to the center of class-four ambiguous islands, etc. Note that since the 
direction of rotation in Eq. reverses, the rays turn in an alternating manner. 

The specification of the symbol boundary is complete once we define the codes for points 
that fall on the boundary. This is especially important for c > 2, as there can be more than 
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Figure 19: Epicycle picture for the class-three orbit ^ : | : |. The boundary of W- (blue) is closed on the solid 
rays and open on the dashed rays. 

two ambiguous islands. For example in Fig. I19L since (72 = 4 is even, there are four ambiguous 
class-two islands, corresponding tot = 4= l:l, 7=1:2, 8 = —1 : —2 and 11 = —1 : — 1 — and 
the rays defining the boundary of W- intersect all of these islands. By contrast, when q2 is odd, 
there are only two ambiguous class- two islands, t = 1 : 1 and —1 : —1. 

The openness/closedness state of the boundary of the region W- is most easily described 
recursively. Recall that at class one, the boundary of W- is open. Given the boundary at class 
(c — 1), the boundary at class c is obtained by reversing the state of the ray at class c — 1 and 
appending class-c rays with the same state. Thus for class-two, the previously open class-one 
ray becomes closed, and the new class-two ray is also closed. For class three, as shown in Fig.lTol 
the class-two ray and the new class-three ray are now open. Applying this rule to the first few 
classes generates Table 0] 



class 


ri 


r2 


r3 


r4 


^5 


1 












2 
3 

4 












5 













Table 4: The class-c boundary of W- is open (dashed) or closed (solid) along the rays, rc, 
defining the ambiguous islands. 

We have checked that this rule works for a number of class-three orbits, and a few orbits 
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of class four. However, the high periods and multipHcity of secondary bifurcations involved 
prohibit systematic studies of high class orbits. In particular for the class-two orbit i : 1, this 

2' 3' 4' 5' 



rule generates the correct codes for the elliptic and hyperbolic orbits with ^3 — ^ ^ ^ ^ 



and |. We have also checked the case \ '■ \ ■ and a number of cases with odd q2 including 
I : ^ : |. The correct code is also obtained for ^ : ^ : ^ : i. The class-c codes also work for 
the period doubling case where it is equivalent to the doubling substitution rule, recall Eq. (|25() . 
For example, the codes for the period-doubling sequence through class four are 



1 . 1 

2 • 2 
1.1.1 
2 • 2 • 2 
1.1.1.1 
2 • 2 • 2 ■ 2 



(- + — r 

(- + — + -+r (31) 
(- + — + - + - + — + — )°° 



A simple rule that also generates this sequence is to double the previous entry and flip the last 
symbol [10] . 

Though we have not systematically explored beyond class-two, these few studies give us 
confidence that the class-c rule works in general. 



6 Conclusion 

We have shown how to systematically generate symbolic codes for rotational orbits and "islands- 
around- islands" orbits for the area-preserving quadratic map. The construction of the codes in 
terms of wedges with opening angles determined by the frequencies uJc of the orbits explains 
numerical observations of the wedge-shaped symbol boundaries for elliptic islands. The wedge 
boundaries are constructed from segments of the symmetry lines of the mapping. 

It is for precisely this case — when there are elliptic orbits — that more traditional, symbolic 
partitions fail. It is interesting to note that Christiansen and Politi constructed symbol bound- 
aries in islands of the standard map that have precisely the same structure as ours for class-one 
rotational orbits 0T]. However, for this case the identification of symmetry lines is complicated 
by the fact that the standard map has two distinct sets of reversors. 

Open questions that we hope will be investigated in the future include: 

• How does the symbol boundary evolve from the line x = when there is a horseshoe, to 
the complex set of epicycle-generated wedges when there are elliptic islands? Note that 
when an orbit first becomes stable, it typically does so by an inverse period-doubling, so 
its wedge opening angle will be zero. 

• Does the epicycle picture for rotational codes also apply to other reversible, area-preserving 
maps with anti-integrable limits? Examples include the standard mapping, and polynomial 
automorphisms . 

• How do the wedge boundaries connect across the island separatrices? As was first observed 
in j2j, the symmetry lines seem to naturally connect with the symbol boundary in the 
chaotic region of phase space that corresponds to primary tangencies. They remark: "We 
have no explanation for this nice phenomenon." Neither do we. 

• What are the codes for rotational orbits generated by bifurcations of the elliptic orbits 
created in secondary pitchfork and twistless bifurcations? 

• How are the symbol boundaries organized in maps that are not reversible? Perhaps the 
simplest example corresponds to some quartic polynomial automorphisms |43j ? 
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A Codes for Twist Maps 



Our purpose in this appendix is to relate the rotational codes described in this paper to those used 
in other papers. There are three commonly used codes for rotational orbits of twist maps, the 
linear, velocity and acceleration codes, which we denote a, b, and c respectively P^IHIHSIUS] . 

Aubry-Mather theory implies that recurrent minimizing orbits for twist maps have a particu- 
larly simple symbolic coding. Orbits that are nondegenerate minima of the action are hyperbolic. 
Thus the codes we discuss here will correspond to the hyperbolic codes in ^ 

An area-preserving map on the cylinder / : x M 

{x',y') = f{x,y) (32) 

is a monotone twist map if ^ > 0. Note that this hypothesis does not apply to the Henon 
map, even locally about (— )°°, because it has a twistless bifurcation at A: = ^ leading to a twist 
reversal j31j . Nevertheless, we will see that the c-codes for minimizing orbits of twist maps do 
correspond to those of the rotational Henon map with the proper translation. 



A.l Linear Code 

One way of coding rotational orbits of maps with an angle variable is to count the number of 
complete rotations they make. To do this, we lift the angle variable to the line, and define the 
linear code, 

at = \xt] ■ 

Thus at is the number of complete rotations at time tM 

The linear code of a minimizing orbit of a twist map is determined by its rotation number. 

Theorem 11 (Aubry [37j). Suppose F : ^ is a lift of the twist map Eq. iT^) . and 

{xt,t G Z} is the configuration of a recurrent, minimizing orbit. Then there exists a rotation 
number uj and phase a such that the linear code of xt is the same as that for the rigid rotation 
on the circle 6t = ojt + 9q, i.e., 

at = \ujt + a] (33) 

Here we take Eq. H33p as the definition of the linear code for a rotational orbit and consider 
some simple structures that arise. 



A. 2 Velocity Code 

An alternative coding for rotational orbits is the "velocity" code, defined to be the first difference 
of the a-code: 

bt = at+i - at (34) 

Note that knowledge of b gives a up to an initial condition, ao, which corresponds to the choice 
of interval for the lifted angle xq. 

While, the 6-code is well defined for any lo, if we restrict to < < 1 then we can compute 
it by dividing the circumference-one circle into two sectors 

Si(o.)EE(-a;,0] , 

^Bo(c^) = (0,1-cu] . 

^^The choice of the ceihng function here is arbitrary, as are the choices of forward and backward finite differences 
in the next two sections. Ultimately, these only affect the canonical ordering for the codes, but the ordering is 
important for concatenation rules. 
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Then we have 

bt = i if {iot + a} e Bi{oj) , (36) 

as shown in Fig. QUI 

When u is irrational, the code changes as a varies. We define the "canonical" ordering as 
that obtained by setting a = 0. For example the canonical 6-code for 7"^ is 

b = (. . . 0100.10100101001001010010100100101001001010010 . . .) . 

The fact that there are uncountably many hyperbolic codes is consistent with the fact that there 
are uncountably many orbits on a particular invariant circle or cantorus. Note that the velocity 
code always consists of blocks of the form (10™) and (10™""^) where m = \_uj~^\ . As we will see 
below, it is no accident that these building blocks are the codes for the Farey parents ^^^q:^ and 
— of each frequency in the interval that they bound. 




Figure 20: Construction of the velocity code for lj = 7 ^ showing the sectors Bo and Bi (shaded). The orbit 
shown has a = -0.1 and its 6-code is (. . . 0100.1010010010 . . .). 

When to = p/q is rational the velocity code is independent of choice of initial phase. 

Lemma 12. The b-code for a rational rotation number uo = p/q is independent of the choice of 
phase a up to cyclic permutations. 

Proof. The Sj defined in Eq. (|15j) partition the circle into q sectors of width ^. Note that 

-Bo = Uj=o ^"^j' ^^'^ -^1 ~ Uj=g-p'S'j so the 6-code is determined by these sectors. Since for 
any a the points 0t = {cot + a} are spaced uniformly on the circle at a distance | from each 
other, there is precisely one 9j in each Sj. Thus the number of and 1 symbols is independent 
of a. Moreover if 9t £ Sj, then 9t+i G Sj+p, so that the cyclic order of symbols is independent 
of a. □ 

Note in particular, that all values —1/q < a < give the canonical velocity code. Other 
values of a, give a cyclic permutation of this "canonical" ordering for the code. 

Corollary 13. For uo = p/q the b-code has p I's and q — p O's. 
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The velocity codes for various frequencies can be easily constructed using a Farey tree pro- 
cedure |231 0] ■ A Farey tree is a binary tree that generates all numbers in an interval between 
two "neighboring" rationals [17|. A pair of rationals | < are neighbors if 

mq — pn = 1 

The Farey tree is recursively constructed beginning with a base defined by a neighboring pair, 
and recursively applying the Farey sum operation 

V m V + m 
-® — = — — , 
q n q + n 

to each neighboring pair. Note that | < | © ^ < and that the daughter is a neighbor to 
each of its parents. 

For example, the two rationals j and j define the base for a Farey tree that includes all 
numbers in the unit interval, see Fig. 1211 The root of the tree, or level zero, is the daughter 
rational defined by the Farey sum of the base rationals. In this case the root is ^. Each 
subsequent level of the tree consists of the rationals that are the Farey sum of each number on 
the previous level with its two neighbors at earlier levels, thus there are 2^ rationals at level j. 

Every irrational to £ [0, 1] is uniquely determined by a path defined by an infinite sequence of 
left, L, and right, R, transitions beginning at the root of the Farey tree. Every rational uj £ [0, 1] 
is uniquely determined by a finite path. For example, the path for 5/13 is LRLR. There are also 
paths in the Farey tree that do not correspond to real numbers: for example LRLLLLLLL . . . 
limits to ^, but is better thought of as Dynamically this sequence corresponds to a ho- 
moclinic orbit. The paths of the parents of any rational are easy to obtain by appropriately 
truncating its path; for example, given a path LRLLL = then the direct parent is obtained 
by simply removing the last symbol, LRLL = The other parent is obtained by removing all 
of the final repeated symbols (L in this case) in the sequence and one more symbol, yielding the 
parent L = ^. 




1/4 2/5 3/5 3/4 

(1000) (10100) (11010) (1110) 




1/5 2/7 3/8 3/7 4/7 5/8 5/7 4/5 

(10000) (1001000) (10100100) (1010100) (1101010) (11011010) (1110110) (11110) 

Figure 21: Farey tree for the base 2 and i for levels zero to three, and the corresponding fe-codes. 
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Lemma 14. The canonical velocity code for the Farey daughter of neighbors ^ < ^ is the 
concatenation of the codes of the parents: 

b -©— =b— b'"^ 



q n J \ n 

Proof. Let lv = denote the daughter, and 9t = ojt. Then ht{uj) = z if {9t} G Bi{oj). We first 
show that the first n symbols are given by the code for — , i.e. that {— i} G Bi{—) imphes that 
{9t} G Bi{uj) for < t < n. 

Suppose first that bt{^) = 1, then there is an integer j such that 



mm , , 

j--<-t<j ■ 37 
n n 

Since a; < the right inequahty imphes ujt < j; providing t > 0. This is the first half of what 
we desired to show. Since t is an integer, the left inequality in Eq. (|S7|) implies that 

m , , 1 

- t+1 > J+- . 

n n 

Combining this with the relation 

p + m m 1 



q + n n n{q + n) ' 
implies that 

1 t + 1 

uj{t + l)>3 + ' 

n n[q + n) 

providing t < q + n — 1. Thus we have shown bt{^) = 1 =^=- bt{uj) = 1, forO<f<g + n — 1 
(this is more than we needed to prove). 

To check that the symbols agree, we must show that when there exists an integer j such 
that 

m m 

j <-t<j + l 

n n 

then {9t} £ Bq{uj). A calculation similar to the previous one shows that this is true when 
< t < g + n. 

So finally we have shown that the symbol sequence b(u;) is given by b(^) repeated, except 
for the last symbol. 

Similar calculations show that all symbols but the first in bt{uj) agree with repeated 
from the end. □ 

It is easy to see that Lemma El is consistent with the codes in the Farey tree in Fig. 1211 For 
example since 6(1/3) = (100) and b(2/5) = (10100), then b(3/8) = b(2/5)b(l/3) = (10100100), 
and since ^ = | |, then b(5/13) = b(2/5)b(3/8) = (1010010100100). 

Since all numbers can be constructed by Farey paths, Lemma extends for irrational 
numbers as well: 

Corollary 15. For any uj whose Farey path includes a parent rational ^ > cu, the first n symbols 
ofb{uj) are those ofb{^). 

Thus, for example the first 34 symbols in the code for 7"^ are given by those of its upper 
Farey neighbor, whose code can be constructed by concatenation: 

6(13/34) = 6(5/13)6(8/21) = 6(5/13)6(5/13)6(3/8) 
= 6(2/5)6(3/8)6(2/5)6(3/8)6(3/8) = . . . 
= (1010010100100101001010010010100100) . 
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A. 3 Acceleration Codes 



The acceleration code is defined to be the magnitude of the first difference of the velocity code 

ct = \bt - bt-i\ . (38) 

This is well defined only for orbits with < uj < ^. It could also be called a "same-different" or 
"exclusive-or" code, since ct = if the velocities are the same and 1 if the velocity changes. The 
point is that for w < i, there is a one-to-one correspondence between allowed b and acceleration 
codes because the 6-code can never have two or more consecutive 1 symbols, thus we obtain 
Q = only for the case of a double in the b code. 

The acceleration code can be obtained geometrically by defining the sectors Ci{uj) = {—u;,lj] 
and Ci{uj) = (a;, 1 — w], so that ct = i when + a} ^ Ci{u;). This follows because we obtain 
Q = 1 whenever 6^ = 1 or bt-i = 1. The canonical 6-code can be reconstructed from the 
canonical acceleration code by using the initial condition bo = 1. 

The geometrical construction shows that the symbol 1 always appears doubled in the accel- 
eration code. Moreover, if m = [w~^J , then the acceleration code consists of blocks of the form 
110"^~^ and 110™~^. The acceleration code inherits properties of the 6-code. 

Corollary 16. The acceleration code for a rational uj is independent of the phase a. 

Corollary 17. The canonical acceleration code for a Farey- daughter is the concatenation of the 
codes for its parents. 

Proof. This follows from Lemma ITU and the fact that the first symbol in the b code is always 1 
and the last is always 0. Thus concatenation does not disturb the calculation of the acceleration 
code symbols. □ 
The interiors of the sectors Ci{uj) are identical to the interiors of the wedges Wi{uj) that 
define the s-code in ^ Recall that hyperbolic codes were those in the interior of the Wi. Thus 
if we translate and 1 into -|- and — appropriately, the acceleration and s-codes are identical 
for hyperbolic orbits. 

Lemma 18. // we set St = sgn(2ct — 1), then the acceleration code for an orbit of rotation 
number lo becomes the hyperbolic s-code. 

For example, since the ^ hyperbolic orbit has acceleration code (110)°°, the correspond- 
ing hyperbolic orbit for the Henon map has s-code ( \-)°°. The elliptic | orbit has code 

( — h +)°°, obtained as usual by flipping the second symbol in the canonically ordered hyper- 
bolic code; however, the elliptic codes do not arise directly from the twist map acceleration 
code. 



B Substitution Rules 

An alternative code for rotational orbits is given by the Farey-substitution rule. This code has 
been used for monotone circle maps [3], and turns out to be identical to the 6-code. 

A substitution rule acts on a symbol sequence s by replacing each symbol with a new se- 
quence. Supposing that s G {L,R}°°, we define two substitution operators for left and right 
transitions: 

Fl{L) = L, Fl{R) = RL 
Fr{L)=RL, Fr{R)=R. (39) 
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Then the symbol sequence for a number whose Farey path begins at a root with symbol sequence 
s, is determined by applying the substitution operators in the same sequence. Thus, since 2/7 
has Farey path LLR starting at 1/2 whose symbol sequence is (RL), the b code for 2/7 is 

FlFlFr{RL) = FlFl{RRL) = Fl{RLRLL) = RLLRLLL . 

Note that this operation is associative, i.e. the above could also be read as 

FlFlFr{RL) = FlFlFr{R)FlFlFr{L) = FlFl{R)FlFl{RL) = FlFl{R)FlFl{R)L = ... 

Lemma 19. The canonical hyperbolic b-code for a frequency uj is given by applying the substi- 
tution rule for its Farey path with the translation L = and R= 1. 

Proof. As shown in Lemma IT^ the canonical 5-code for a Farey daughter is obtained by con- 
catenation of the codes its parents. We will show that this concatenation rule is also valid for 
the Farey-substitution rule. Combining this with the fact that the code for ^ is RL = 10, gives 
the result. 

As discussed in ADDendix lA.2l the Farey paths for the two parents are obtained by truncating 
the path of the daughter by removing the last symbol, and all of the repeated last symbols plus 
one more, respectively. Thus there are two possible cases. First if the path for the daughter is 
WRL^ for an arbitrary sequence W and n > 1, then we must show the concatenation rule 

WRL'^ = W® WRL""-^ . 

Translating this into substitutions applied to the root (RL), and using the associative property, 
we obtain 

FwFrF£{RL) = FwFrF2-\RLL) = FwFrF^-\RL)WFrFI~\L) 
= FwFrF2-\RL)Fw{RL) = iW(B WRL""-^) , 

which is what we wanted to show because the © concatenation writes down the code of the right 
neighbor first. A similar calculation holds in the other case, where the daughter is WLR'^, and 
we can show 

WLR^ = WFLR"-^ © W . 

□ 

There are other substitution rules that apply to period doubling, and more generally n- 
tupling sequences 0. For example, the doubling substitution rule is 

LR and R^ LL 

Applying this sequence to the sequence of w gives the sequence for 1/2 : uj. Applying this 
recursively gives the lo bifurcations of any orbit from the period-doubling sequence. 

C Symmetries and Codes of the Horseshoe 

In this section we give simple geometric construction of the relation between the symbol plane 
and an idealized phase space for the horseshoe map. This helps to visualize the relation between 
the symmetry lines and the symmetries present in the symbolic codes. Listead of the standard 
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smooth horseshoe map, we will use a discontinuous, area-preserving horseshoe map that is 
defined everywhere on the square [—1,1] except for the y axis: 



H{x,y) 



(-2x-l,-i(l + y)) , rc<0 
(2x-l,i(l+y)) , x>0 



see Fig. 1221 The region x < corresponds to the symbol .— , and x > to +. Thus the fixed 
points of H are (+)°^ = (1, 1), and (-)°° = (-i, -i). 




Figure 22: Idealized, discontinuous horseshoe map. 



This horseshoe map is reversible, with symmetry S{x, y) ■ 
as H = RS, where 

' (-2y-l,-l(l + x)) 
(2y-l,i(l+x)) , 

{y = -i(l+x) : 



R = HS 



Thus the fixed set of R is given by 7^ = 
Alternatively, we can factor H = ST as 



= (y, x). As such, it can be factored 
y <0 

y < 0}U{y = i(l + x) : 2/ > 0}. 



T = SH 



(-i(l + y),-2x-l) , x<0 
(i(l + y),2x-l) , x>0 



with a fixed set T = {y = — 2x — 1 : x < 0} U {y = 2x — 1 : x > 0}. As usual, we divide the 
symmetry lines into rays at the "elliptic" fixed point (— )°°, with the subscript i denoting the 
rays that lead to the hyperbolic point and o those that lead away, see Fig. 1231 

We also show in Fig. [231 the period-two point (— .-|-)°° = (i,— |) and its image (-|-.— )°° = 
(— |, i). Note that these points lie on % and 7^, respectively, in accord with Table|21since HTq 



e Si, (+ + 

3 orb 
f ) G So. 



o 
oo 



IZo- Similarly, the elliptic period 3 orbit has points ( — h .+ 

("Ij i) £ '?05 and (H — .+)°° = (|, —5)- The hyperbolic period 3 orbit has points ( .-|-)'^ 

(f , -|) + - -r = (-7' I)' and (+ -.-)- = (-f , 

We can divide the symmetry lines into rays at the (— )°° fixed point. Note that points with 
X larger than that for the elliptic fixed point will be on the ray that goes to the hyperbolic 
point (labeled i), and those with x smaller will be on the ray o. A standard result gives us the 
ordering. 

Lemma 20 (Ordering (c.f. [3", §2.3.2])). Suppose two symbol sequences s = .sqSiS2--- and 

s' agree for their first j symbols, but that sj 7^ s'^. Then the corresponding points x and x' on 
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Figure 23: Symmetry lines and periodic orbits of the discontinuous horseshoe. The two fixed points are shown 
as a circle and x, the period-two points as squares, and the points on the two period-three orbits as triangles. 

the horseshoe are ordered as x < x' if the parity, Eq. of sqSi ■ ■ ■ sj is odd (and hence that of 
sqSi ■ ■ ■ s'j is even). 

Thus sequences are ordered the same as their parities. This imphes that symbol sequences for 
the symmetry rays are distinguished by "odd" or "even" blocks of — symbols after the binary 
point, see Eq. (j2T|) . 
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